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1.0  INTRODUCTION 


This  research  project  is  directed  towards  examining  the  impact  of  gas 
injections  from  large  space  structures  operating  in  low  earth  orbit.  Our 
particular  emphasis  is  on  phenomena  which  may  lead  to  enhanced  plasmas  about 
the  structure. 

This  first  annual  report  on  the  project  describes  several  related  efforts 
aimed  to  develop  an  understanding  of  the  relevant  plasma/flow  phenomena  with 
the  ultimate  goal  of  designing  a  Shuttle-borne  experiment  to  study  plasma 
enhancement . 

One  of  the  primary  Shuttle  diagnostics  for  studying  plasma  properties  is  the 
Langmuir  probe.  Section  2  provides  a  review  of  probe  theory  with  a 
re-analysis  of  some  of  the  measurements  performed  on  earlier  Shuttle  flights. 

A  relatively  detailed  description  of  our  modeling  on  critical  ionization 
velocity  (CIV)  effects  is  provided  in  Section  3.  The  critical  ionization 
velocity  effect  originally  proposed  by  Alfven  (195A)  may  be  an  important 
source  of  plasma  enhancement  around  large  space  structures  such  as  the  Shuttle 
orbiter  or  the  future  space  station.  Papadopoulos  (1984)  proposed  that  the 
glow  observed  near  the  surfaces  of  the  Shuttle  was  the  result  of  electrons 
heated  by  a  plasma  instability  produced  by  reflection  of  some  of  the 
ionospheric  ions  when  the  Shuttle  moves  through  the  ionosphere.  Indeed,  the 
ionic  environment  about  the  orbiter  is  observed  to  have  periods  where  the  ion 
density  increases  dramatically  (Siskind,  et  al.,  1984).  The  electron  energy 
distribution  often  is  not  a  simple  Maxwellian  distribution  (Raitt,  et  al., 
1984;  Murphy,  et  al.,  1986;  Raitt,  et  al.,  1987).  In  addition  orbiter 
operations,  such  as  thruster  firings,  flash  evaporator  releases,  and  water 
dumps  are  observed  to  enhance  electrostatic  noise  and  to  increase  the  plasma 
density  fluctuations  (Pickett,  et  al.,  1985).  Thus  there  is  evidence  that 
plasma  Instabilities,  and  possibly  CIV  processes,  are  produced  as  the  result 
of  orbiter  operations.  If  CIV  processes  are  occurring,  it  is  not  clear  if 
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these  result  from  reflection  of  streaming  ions  (Papadopoulos ,  1984), 
ionization  of  metastable  states  (McNeil,  Lai,  and  Murad,  1987),  or  ionization 
of  dimers  or  trimers  formed  during  the  free  expansion  into  space  (Murad,  et 
al.,  1986).  Section  3  describes  a  model  for  the  two  stream  instability  which 
is  used  to  gauge  the  importance  of  CIV  around  Shuttle. 

The  primary  objective  of  this  program  is  to  conceptually  design  a  well 
diagnosed  injection  experiment  to  be  flown  on  space  Shuttle  for  the  purpose  of 
studying  the  critical  ionization  velocity  process.  This  experiment  will 
concomitantly  provide  reliable  data  on  the  "ambient"  plasma  fields  around  the 
orbiter.  The  experiment  would  employ  infrared  active  species  and  be  deployed 
in  conjunction  with  the  IBSS  experiment  in  order  to  take  advantage  of  the 
outstanding  ultra-violet  and  infrared  radiation  diagnostics  available  in  that 
payload.  The  CIV  experiment  would  also  utilize  radiation  and  plasma 
diagnostics  placed  within  the  Shuttle  bay. 

Specifically  we  propose  to  inject  gases  form  the  Shuttle  bay  in  low  earth 
orbit  under  conditions  which  should  lead  to  enhanced  plasma  due  to  the  CIV 
phenomena.  The  key  is  in  choosing  the  right  species  to  inject.  Gases  cool 
when  they  expand  into  a  near  vacuum,  as  the  random  motion  in  thermal  energy  is 
converted  into  directed  motion.  The  cooling  gas  forms  dimers,  trimers,  and 
possibly  higher  clusters.  A  dimer  (or  a  higher  cluster)  has  two  advantages 
over  the  corresponding  monomer  in  reducing  its  critical  velocity  (v^.)  value: 
it  has  a  slightly  lower  ionization  potential  and  it  has  twice  (or  n  times)  the 
mass.  Thus  releasing  the  gas  under  conditions  where  the  free  expansion  favors 
production  of  dimers  or  trimers  should  enhance  the  prospects  for  CIV  processes 
(Murad,  et  al.,  1986).  Our  experimental  gas  release  will  form  such  dimers  and 
higher  clusters,  and  their  smaller  Vj,  values  will  enhance  the  probability  that 
the  CIV  process  will  produce  enhanced  ionization. 

At  present  we  intend  to  inject  four  distinct  gases.  First  a  null  injection, 
utilizing  a  rare  gas  such  as  neon.  This  species  does  not  have  a  favorable 
ionization  potential  for  fostering  CIV  at  low  earth  orbital  velocities  of  8 
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km/s  and  thus  any  observed  plasma  fluctuations  may  be  ascribed  to 
perturbations  due  to  the  injected  flow.  Two  other  gases  to  be  injected  are  NO 
and  002-  These  species  have  both  the  merit  of  being  typical  rocket  exhaust 
gases,  and  having  dimers  with  ionization  potentials  corresponding  to 
equivalent  critical  ionization  velocities  below  orbital  velocity.  Indeed  for 
ram  direction  injection  the  monomers  themselves  may  ionize  through  CIV.  These 
polar  molecules  and  their  ions  exhibit  infrared  bands  which  are  within  the 
spectral  range  of  the  IBSS  payload.  Furthermore  the  bands  of  the  neutral  and 
ionic  species  may  be  resolved  from  each  other.  Lastly  in  the  case  of  C02''', 
the  dissociative  recombination  process  produces  electronically  excited  CO, 
whose  radiation  can  be  monitored  by  the  AIS  instrument.  Thus  any  radiation 
resulting  from  the  CIV  effect  directly  caused  by  the  injection  of  these  two 
species  may  be  readily  monitored  by  the  IBSS  payload.  The  last  species  to  be 
injected  will  be  xenon  which  should  be  readily  ionized  by  the  CIV  process 
(Mobius,  et  al.,  1979;  Axnas,  1980). 

The  experiment  design  is  still  in  the  preliminary  stages.  A  Phase  0/1 
Accident  Risk  Assessment  report  has  been  filed  and  is  included  as  an  appendix 
to  this  volume.  Furthermore  a  Preliminary  Design  Review  has  just  been 
completed.  The  viewgraphs  presented  at  this  review  are  provided  in  a  separate 
volume. 
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2.0  MEASUaEMENT  OF  ELECTRON  TEMPERATURE  AND  DENSITY 


The  measurement  of  electron  density  (n^)  and  temperature  (T)  in  the  Shuttle 
environment  can  readily  be  performed  by  use  of  a  Langmuir  probe.  The  theory 
of  Langmuir  probes  under  conditions  where  the  particle  mean  free  path  is  large 
as  compared  to  the  probe  and  sheath  dimensions  and  where  the  effect  of 
magnetic  field  on  charged  particle  orbits  is  negligible  is  well  understood 
(Langmuir,  1961;  Chen,  1963;  Schott,  1968).  We  show  in  Fig.  1  the  plasma 
sheath  thickness  as  a  function  of  ng.  Since  a  plasma  is  a  good  conductor,  the 
electric  field  will  be  zero  inside  the  plasma  and  application  of  a  voltage 
onto  a  probe  in  contact  with  the  plasma  will  result  in  a  current  being  drawn 
from  the  edge  of  the  sheath.  The  electron  gyroradius  in  a  magnetic  field  (B) 
of  0.3  Gauss  is  given  by: 

R  =  11.2  ^If  (cm)  (1) 

where  T  is  in  eV.  R  varies  inversely  with  B.  The  ion  gyroradius  will  be 
larger  by  the  square  root  of  the  ion  to  electron  mass  ratio.  One  sees  from 
Eq.  (1)  that  if  the  probe  dimensions  are  of  the  order  of  3  cm  or  less,  then 
electrons  will  have  straight  orbits  and  magnetic  field  effects  can  be 
neglected.  At  low  electron  densities  (ng  <  10^  cm“3)  the  probe  when  biased 
positively  with  respect  to  the  plasma,  will  draw  electrons  from  a  distance  so 
large  that  the  magnetic  field  is  expected  to  play  a  role  in  the  conduction 
independent  of  probe  size. 

We  show  in  Fig.  2  a  typical  current  versus  voltage  characteristic.  In  region 
III  the  probe  is  biased  negatively  with  respect  to  the  plasma  and  is  attract¬ 
ing  ions.  In  region  I  the  probe  is  biased  positively  and  is  attracting  elec¬ 
trons.  When  the  probe  is  at  the  sheath  potential  Vg,  i.e.,  when  there  is  no 
voltage  difference  between  the  probe  and  the  plasma,  the  current  to  the  probe 
is: 
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where  Ap  is  the  area  of  the  probe,  and  the  mean  velocity  (8  of 

the  ions  (electrons).  In  the  last  step  we  used  the  fact  that  the  ion  current 
could  be  neglected  since  the  electron  thermal  velocity  is  so  much  larger  than 
the  ion  velocity.  The  floating  potential  Vf  is  the  potential  at  which  the  ion 
and  electron  currents  cancel  and  where  the  probe  draws  no  net  current,  the 
transition  region  TI  is  a  region  where  the  probe  is  biased  negatively  with 
respect  to  the  plasma  but  not  so  much  as  to  completely  repel  the  electrons. 

It  is  a  region  where  the  current  is  still  dominated  by  flow  of  electrons  to 
the  probe. 


The  maximum  current  that  a  probe  can  collect  is: 
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where  np  is  the  plasma  density  (np  =  n^)  and  Ag^  the  area  encompassed  by  the 
edge  of  the  sheath.  The  slope  of  the  I  versus  V  curve  in  region  II  will  yield 
the  electron  temperature,  since,  when  electrons  are  repelled,  it  can  be  shown 
that  Eq.  (1)  still  holds  with  Ag^  being  equal  to  the  probe  area  Ap  and  n 
varying  as 


n  =  ng  exp(-  e  |v  |/kT) 

where  V  is  the  (negative)  potential  with  respect  to  the  plasma.  We  therefore 
have 


d  In  I  e 

=  ■  RT 


(4) 


Once  the  electron  temperature  is  known,  the  electron  density  can  be  derived 
from  the  saturation  current  which  varies  proportionately  with  np'Tf.  One 
difficulty  in  determining  np  exactly  is  that  the  sheath  area  will  change  with 
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applied  voltage.  When  the  sheath  is  very  thin  as  compared  to  the  probe 
dimensions  then  Ag^  is  practically  equal  to  the  probe  area  Ap,  and  the 
saturation  current  appears  to  be  practically  independent  of  probe  voltage. 

When  the  probe  is  attracting  electrons  one  can  calculate  the  current  by 
following  the  electron  trajectories,  starting  from  the  edge  of  the  sheath. 

The  geometry  is  shown  in  Fig.  3,  for  either  a  spherical  probe,  or  an 
infinitely  long  cylindrical  probe.  A  particle  starting  at  the  edge  of  the 
sheath  at  point  M  will  hit  the  probe  (at  point  N)  if  the  angle  0  shown  in  the 
figure  is  less  than  some  critical  ang’o  0^  which  depends  on  the  initial  radial 
velocity  u.  By  integrating  over  all  u  and  all  angles  0  <  Oq(u)  one  obtains, 
for  a  spherical  probe,  the  current: 

I  =  Ap  jr  F  (5) 

where  Ap  =  Aria^,  =  (ngC/4)  and 

2 

F  =  ^  (1  -  exp(-  ♦])  +  exp(-  (6) 

a 

with 


If  the  sheath  thickness  is  very  small  as  compared  to  a,  then,  when  eV  >  kT, 

♦  »  1,  F  =  s^/a^  =  1,  and  the  current  is  given  by  Eq.  (3)  with  Ag^  =  Ap.  If 
T  is  known,  then  it  is  an  easy  matter  to  derive  n  from  the  current.  If  the 
sheath  thickness  is  small  as  compared  to  a,  we  can  use  the  Langmuir-Tonks 
equation  for  space  charge  limited  flow  in  planar  geometry  to  calculate  the 
sheath  thickness  as  a  function  of  probe  voltage  with  respect  to  the  plasma 
(Coblne,  1968). 
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Figure  3.  Electron  trajectory  in  the  sheath. 


/>)  ,,3/2  ,  A 

I  =  A  - V -  =  2.33  X  10"^  -  P— , 

^  9n(s-a)  (s-a) 


,3/2 


(amps) 


(8) 


where  in  the  last  step  Ap  is  in  cm^,  v  in  volts  and  (s-a)  in  cm.  Both  s  and  n 
can  be  calculated  from  Eqs.  (3)  and  (8)  if  T  is  known  since  Agjj/Ap  =  s^/a^. 


When  the  sheath  thickness  is  comparable  to  or  larger  than  the  probe  radius,  it 
is  not  straightforward  to  derive  n  from  the  saturation  current,  and  one  must 
solve  Poisson's  Equation  taking  into  account  particle  trajectories  with  plasma 
density  as  a  parameter.  Only  a  detailed  comparison  of  the  theoretical  I 
versus  V  characteristics  with  the  measured  one  will  yield  n,  once  T  is  known. 


2.1  ANALYSIS  OF  SHUTTLE  PROBE  DATA 


Murphy,  et  al.  (1986)  used  a  spherical  probe  of  radius  a  =  3  cm  to  measure  n 
and  T  in  the  Shuttle  environment.  They  present  two  characteristics 
corresponding  to  n  =  2  x  10^  cm"^  and  n  =  lO'^  cm“3.  Their  data  is  shown  in 
Figs.  4  and  5.  They  calculated  Tg  by  fitting  the  slope  given  by  Eq.  (4)  to 
region  II  of  the  characteristic,  after  correcting  the  characteristic  for  the 
ion  current.  Because  of  the  large  velocity  of  the  Shuttle  with  respect  to  the 
ionospheric  plasma,  the  ion  current  will  be  Ij  =  njWiia^,  where  W  (=  8  km/s)  is 
the  Shuttle  velocity  with  respect  to  the  plasma.  This  is  an  order  of 
magnitude  larger  current  than  one  would  draw  if  the  probe  had  no  motion  with 
respect  to  the  background  plasma.  Based  on  the  electron  temperature  and 
density  that  they  derive  from  their  data,  we  calculate  the  sheath  thickness  to 
be  3  and  0.2  cm  for  Og  =  10^  and  2  x  10^  cm~3,  respectively.  The  sheath 
thickness  is  therefore  of  the  order  of  the  probe  radius  at  the  lower  density. 
We  have  used  Eqs.  (5)  to  (7)  to  calculate  the  electron  current  as  a  function 

10,000 
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1 

Figure  4.  Current  versus  voltage  n  =  2.6  x  10^  cm“3. 
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PROBE  VOLTAGE  (Volts) 

A-5326 

Figure  5.  Current  versus  voltage  n  =  10^  cm“3. 

of  V  for  V  positive.  The  potential  of  the  probe  is  measured  with  respect  to 
the  Shuttle  body,  however,  and  the  potential  with  respect  to  the  plasma  at  the 
location  of  the  probe  is  unknown.  We  translated  the  characteristic  so  that 
the  calculated  current  would  fall  on  top  of  the  measured  current  in  the 
retarding  region  II.  We  took  s-a  =*  3  cm  when  ng  =  10^  cm~^  and  calculated  s 
from  Eqs.  (3)  and  (8)  when  ng  =  2  x  10^  cm~^.  We  see  that  we  obtain 
reasonable  agreement  with  Murphy  et  al.'s  measurements  when  ng  =  2  x  10^  cm~^, 
but  find  a  saturation  current  too  large  by  a  factor  of  three  when  we  set 
ng  =  2  X  10®  cm~^.  Our  results  indicate  that  the  actual  plasma  density  is  a 
factor  of  three  to  four  times  lower  than  the  density  derived  by  Murphy  et  al. 

2.2  EFFECT  OF  SOLAR  RADIATION  ON  MEASURED  CURRENT 

Photoelectrons  are  generated  by  the  solar  XUV  photons  that  are  incident  in  the 
probe  and  will  affect  the  current  when  the  probe  is  biased  negatively  with 
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respect  to  the  plasma.  We  calculate  the  effect  for  a  gold  plated  probe,  such 
as  the  one  used  by  Murphy  et  al.  (1965).  Let  tgCX)  be  the  solar  flux  in 
photons/cm^  s  per  unit  wavelength  Interval.  Let  r|(X)  be  the  quantum 
efficiency  of  the  electron  ejection  process.  Integrating  over  the  probe  and 
over  wavenumber,  we  calculate  the  following  photo-electron  current  (q  = 
electron  charge). 


n/2 

=  -2iiqa^  J  cos  6  d(cos  0)  J  ♦g(^)  h(X)  dX 
0 

Ig  =  nqa^  J  h(X)  dX  . 

We  show  in  Fig.  6  the  photo-electric  yield  of  gold  (Samson,  1967).  It  only 
becomes  appreciable  at  wavelength  X  <  130  nm.  Integrating  the  photon 
efficiency  times  an  average  solar  flux  over  X,  we  find  a  photo  electric 
current  of  -8  x  10^  el/cm2  s  leading  to  Ig  =  3.6  x  lO-^A.  This  current  is  a 
factor  of  four  below  the  detection  threshold  in  Murphy  et  al.s  (1968) 
measurements. 


Figure  6.  Photoelectric  yield  of  Au  as  a  function  of  wavelength. 
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2 . 3  CONCLUSIONS 


Probe  theory  under  the  conditions  of  Murphy  et  al.'s  (1968)  experiments 
reproduces  pretty  well  the  experimental  results  at  high  electron  densities 
(n  =  2  X  10^  cm“3)  but  indicates  that  at  lower  electron  density  (n  =  10^ 
cm~3),  where  the  sheath  thickness  is  comparable  to  the  probe  radius,  they 
overestimated  n  by  a  factor  of  three  to  four.  The  photo-electric  yield  of 
gold  is  shifted  to  short  enough  wavelengths  that  the  current  generated  by 
solar  illumination  is  below  the  measurement  threshold. 
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3.0  SIMULATION  OF  CRITICAL  IONIZATION  VELOCITY  EFFECTS 


The  critical  ionization  velocity  (CIV)  effect  was  first  proposed  by  Alfven 
(1954)  in  the  context  of  his  model  on  the  origin  of  the  solar  system.  The 
proposal  states  that  there  is  a  rapid  ionization  whenever  a  neutral  gas  moves 
across  the  magnetic  field  with  a  velocity  relative  to  a  plasma  exceeding  the 
value 

Vc  =  (2I/rab)0-5  , 

where  I  is  the  ionization  potential  of  the  neutral  gas  and  mjj  is  the  mass  of 
the  ion  formed  upon  ionization.  There  is  now  considerable  evidence  for  CIV 
ionization  at  velocities  near  v^.  from  various  laboratory  experiments 
(Danielsson,  1973;  Himmel  et  al.,  1976).  Recently  there  have  been  experiments 
with  chemical  releases  in  space  (Newell,  1985).  One  of  these  shows  CIV 
ionization  (Haerendel,  1982),  but  others  show  only  a  small  amount  or  do  not 
detect  any  (Wescott,  et  al.,  1986a, b;  Kelley,  et  al.  1986;  Torbert  and  Newell, 
1986).  In  other  releases  the  plasma  changes  could  result  from  processes  other 
than  CIV  (Sasaki,  et  al.,  1985),  or  from  the  ion  beam  released  (Sasaki,  et 
al.,  1986). 

We  have  investigated  CIV  phenomena  for  conditions  that  can  occur  when  a  beam 
of  neutral  gas  molecules  is  released  from  the  Shuttle  in  low  earth  orbit. 

Here  we  are  interested  in  modeling  both  the  molecular  phenomena  and  the  plasma 
effects  that  result  from  the  formation  of  a  beam  of  fast  ions  moving  across  a 
magnetized  plasma.  We  first  present  the  the  linear  dispersion  theory,  then  we 
present  simulation  results  for  a  simplified  plasma.  Next  we  present  more 
linear  dispersion  theory  for  a  warm  plasma.  We  conclude  with  simulations  of 
the  molecular  and  plasma  effects  in  combination. 
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3.1  LINEAR  DISPERSION  THEORY  OF  WAVES  IN  A  "COLD"  PLASMA 


3.1.1  Introduction 

Here  we  are  concerned  with  small  amplitude  wave  imstlons  of  a  cold  Ionized  gas. 
This  Is  one  of  the  few  plasma  situations  for  which  It  Is  possible  to  give  a 
systematic  account.  The  cold  plasma  Is  of  Interest  because  It  can  support  many 
of  the  varied  types  of  wave  motion  which  are  present  at  non-zero  temperatures, 
but  without  any  of  the  mathematical  subtleties  that  are  Inherent  In  a  solution  of 
the  kinetic  equations  for  a  hot  plasma.  A  cold  plasma  means,  In  this  context,  an 
unbounded,  electron-ion  plasma  of  zero  thermal  velocities.  Such  a  system  Is 
actually  a  model,  the  limits  of  the  validity  of  which  can  best  be  Inferred  from 
the  theory  for  a  hot  plasma. 


In  what  follows  we  develop  the  linear  theory  of  waves  In  a  “cold"  plasma.  This 
follows  closely  the  text  of  Montgomery  and  Tldman  (1964).  The  theory  Is 
subsequently  applied  to  a  situation  of  practical  Interest  In  the  design  of  the 
space  shuttle  experiment  to  test  the  critical  ionization  velocity  hypothesis. 


Let  us  assume  that  each  of  the  species  of  charge  has  a  number  density  n^,  a 
velocity  V£,  no  pressure  gradients,  and  obeys  an  equation  of  continuity  and  an 
equation  of  motion: 


3 

n"l 


+  V»(njV^) 


0 


(9) 


3t 


®1  1  -  - 
V  V,  =  —  (E  +  -  V,  X  b) 
1  m^  ^  cl  ' 


(10) 


E  and  B  satisfy  Maxwell's  equations 


7  •  E  =  4n  p 


7  X  E  - 


i  1® 

c  3t 


7  •  B 


7  X  B 


0 

^  -  J.  3| 
c  ^  c  3t 


(11) 

(12) 


(Note  that  we  are  using  Gaussian  units  here.)  The  charge  and  current 
densities  are: 


?  'i-i 


I 

1  ^ 


(13) 
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3.1.2  Perturbation  Theory 


The  Eqs.  (9)  through  (13)  with  appropriate  Initial  and/or  boundary  conditions, 
determine  the  plasma  variables  of  state  nj^  and  V£  as  functions  of  x  and  t. 

They  are  nonlinear  and  very  difficult  to  solve.  A  perturbation-theoretic  solu¬ 
tion  about  a  uniform  equilibrium  Is  possible,  however,  and  useful  for  many 
purposes. 


Let  us  assume  a  unlfrom  magnetic  field  Bq  In  the  zero  order  system  and  set 
-(1)  -  -(1) 

V=u+v  E  =  E 

^  ^  ^  (lA) 


-  -  -(1) 
B  =  B  +  B 
0 


(1) 

n  =  N  +  n 
111 


where  u^,  Bg,  and  Nj  are  Independent  of  x  and  t. 


Let 


J^Ne  =  0  ,  5[Neu  =  0 

11  111 
1  1 


(15) 


and  u  X  B  =0  (16) 

1  0 

l.e.,  there  Is  no  charge  or  current  density  the  In  zero-order  system  and  the 
zero  order  velocities  are  parallel  to  Bg  (no  force). 


Substituting  Eq.  (14)  Into  Eqs.  (8)  through  (13)  and  keeping  linear  terms  only 
gives  the  linear  system  of  equations: 


where 


9  + 

"i  ^ 

N^7  • 

+ 

•  Uj^  =  0 

(17) 

+ 

at  1 

“i  * 

=  !i[E<^> 

m^ 

+  -  vj^^  X  B  +  -  u  X  B^^^  ] 
cl  0  c  1  J 

(18) 

7  X  B^^^ 

4n  - 

(1)  ^ 

1  3  pd) 

,  5  •  •  0 

(19) 

V  xE('> 

1 

a  «  ^  ' 

c 

.  V  . 

(20) 

»  1 

1 

ejN^v 

U), 

ni^>ui) 

(21) 

-  I 

11 

(22) 
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The  linear  system  (17)  through  (22)  has  constant  coefficients,  and  may  be 
solved  using  Fourier  transforms  in  x  and  t.  Hence  we  seek  solutions  for  which 
all  the  perturbed  quantities  have  a  dependence  2,^^^  ^  so  that  (17) 

through  (22)  become  an  algebraic  system  of  linear,  homogeneous  equations. 
Initial  and/or  boundary  conditions  can  be  satisfied  by  linearly  superposing 
these  solutions  into  wave  packets,  l.e.,  Fourier  distributions  In  k  and/or  o). 


Thus,  we  define  the  symmetric  Fourier  transform  pair: 

F(x,t)  =  — /  d\da)  F(k,a))  e^^*^  *  x  -  wt)  ^23a) 

(2it)^ 

F(k,aj)  = — ^  f  d^xdt  F(x,t)  e"^^^  •  x  -  wt)  ^23b) 

(2ir)^ 

Time  differentiation,  taking  the  dlv.  and  taking  the  curl  in  the  direct  space 
Is  equivalent  to  multiplying  the  Fourier  coefficients  by  (-luj),  ik,  and  Ikx 
respectively.  The  relations  between  the  Fourier  coefficients  for  (17) 
through  (22)  are: 


-loanj^^  +  lie  •  k  •  u^  =  0 


(17') 


-luivj^^  +  lu  •  k  vj^^  =  —  +  -  vj^^  x  B  +  -  u  X  ] 

1  1  im^'-  cl  Ocl  (18') 


lie  X  B^^^  =  —  ,  lie  •  B^^^  »  0 

c  c 


ik  X  b(‘>  .  Ik  .  E<‘>  .  4. 

c 


~(1)  r  ^  ~(1) 

P  “I  e.n 
1 


(19*) 


(20') 


(21') 


(22') 
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Equations  (17*)  through  (22*)  are  a  closed  set  for  the  Fourier  coefficients  and 
could  be  used  to  solve  for  them,  whence  by  integration  (inverse  Fourier 
transform),  our  original  system  of  equations  could  be  solved  for  the  plasma 
variables  of  state  (under  the  appropriate  boundary  conditions).  Having  assumed 
a  normal  mode  solution,  the  linear,  homogeneous  system  (17*)  through  (22*)  gives 
us  the  normal  modes  (eigenvectors)  and  normal  frequencies  (eigenvalues). 

In  the  next  section  we  shall  be  interested  in  a  case  for  which  uj  x  Bq  /  0. 

The  derivations  presented  here  apply  provided  we  do  not  consider  {c  =  w  =  0. 

This  is  not  a  serious  limitation  since  we  shall  not  be  interested  in 
transforming  back  to  the  direct  space  (x,t). 

Our  approach  is  to  eliminate  all  perturbed  variables  In  favor  of  the  perturbed 
electric  field.  In  this  form  we  will  obtain  the  dot  product  of  a  dyadic  with 
the  perturbed  electric  field  equal  to  zero.  The  dyadic  may  be  Interpreted  as 
the  dielectric  response  of  the  plasma.  Taking  the  curl  of  (20*)  we  find 


Ik  X  ( Ik  X 


Ik  X 

c 


k  .  v^^^ 

lur^TT  V  rx,  ~(1)  ,  XT  f  “  ^ 

-[rJ'JVi  -""i — “i) 

1  u)  -  k  •  u. 


c 


In  arriving  at  this  result  we  have  used  (19*),  (21*),  and  (17*),  respectively. 
This  equation  can  be  simplified  by  using  the  vector  Identity 

a  X  (b  X  c)  =  (a  •  c)b  -  (a  •  b)c 
The  resulting  simplified  equation  is 


.  2  2,2.  ~(1)  ^  2-x-  ^(i)x  ,  ^  -(1) 

((jj  -  ck)E  +c  k(^k  •  E  J  =  -Airiuj 

r  ~(1) 

=  -4Trlaj  I  e  N  (v,  +  u  — 

1  (oj  -  k  •  u 


. 


(24) 


On  the  right  hand  side  v^  is  still  present.  However,  Eq.  (18*)  (and  Eq.  (22*)) 
are  still  at  our  disposal.  Introducing  the  cyclotron  frequency  vector 
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(25) 


e.B 

5  .-Lo 

1  m^c 


Eq.  (18')  becomes: 


-1(a)  -  ic  •  u^) 


.:i  [E<‘>  .B»>]  - 

m.  '■  c  l  ■' 


1 


In  this  last  equation,  =  -^  it  x  (see  Eq.  (20'))  and  hence  It  remains 

to  solve  for  In  terms  of  E^^\  We  do  this  by  first  recalling  the  Identity: 


^l(^l  * 


n^x  (iij  X 


Next,  we  take  the  dot  and  cross  product  of  Eq.  (26)  with  Oj;  trying  to  take 
advantage  of  Eq.  (27).  From  the  dot  product  of  Eq.  (27)  with  Oi  we  get 


?(1)  ^  1  n 


-1(a)  -  k  •  Uj^) 


•  (u^  X  B^^O)  . 


Adding  this  result  to  the  curl  of  Eq.  (26)  with  flj  we  find 


_  ^(1) 

-Ko)  -  •  u  )  ^  +  - — -  (ii^  •  i  “l  •  ^“l 

^  nj  -Ko)  -  k  .  u^)  ^  ^  ^  0^ 


®i  ^1  ,~(i)  1  -  ~(i)^  ”i^”i  *  ''i 

- s-x^E  +-u  xB  J+ - 5 - 

"i 


Sj  X  (iij  X 
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- 3—2  {!(<*)  -  k  •  u^)  +  i  u^  X  (ic  X  e''^^)] 

-  (o)  -  k  •  Uj^) 


+  jj  X  X  (k  X  E^^^)] - - - 

1  to  1  '•  '  ^  ,  r  ~ 


(to  -  k  •  Uj^) 


•  ii  •  +  —  u,  X  (k  X  E^^^)]  .  (30) 

i  *■  u)  1  '•  '  ^ 

Substitution  of  Eq,  (30)  Into  (24)  gives  an  equation  of  the  form 

R  •  E^^^  =  0  ,  (31) 

♦* 

where  R  Is  a  3  x  3  dyadic  Involving  known  quantltltes.  R  may  also  be  written 
as  a  3  X  3  matrix.  It  Is  well  known  that  a  necessary  condition  for  such  an 
equation  to  have  a  nontrivial  solution  Is  that  the  determinant  of  the  matrix 
vanish.  Letting  |r|  =  0  gives  an  equation  for  the  normal  frequencies 
(0  =  (o(k) . 
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3.2  ION  BEAM  STREAMING  INTO  A  BACKGROUND  ION-ELECTRIC  PLASMA  IN  THE  COLD 
PLASMA  APPROXIMATION 


Let  us  fix  a  coordinate  system  on  the  background  Ion-electron  plasma  and  take 

A 

B  =  (0,0, Bo)  “  BqZ  (32a) 


k  -  (kv,0,k2)  =  kjjX  +  k^Z 


(32b) 


u  =  (V,0,0)  *  VX 


(32c) 


We  want  to  solve  the  normal  mode  problem  for  a»(k).  Specifically,  we  solve 
(see  Eq.  (24)  of  previous  subsection) 


(ui^  -  c^k^)  E  +  c^k(k  •  E)  =  -4iTiujJ 


(33) 


where 


i  -  I  ejN  (I  - - z^-r-) 

i  (u)  -  k  •  u^) 


(34) 


and 

with 


'"i  =  ''id)  ''i(2) 


(35a) 


u  -  n2 

-  (u)  -  k  •  u^) 


{i(oj  -  k  •  Uj^)  E  +  X 


E  -  i 


(u  -  k  •  u^) 


•  E} 
(35b) 


■®i/™i 


0. 


-()-  _  ^-2  {[i(a)  -  k  •  u  )  +  n  X  -  i - 

-  (o)  -  k  •  u^)^  ^  ^  ^ 


(u  -  k  •  Uj^) 


u^  X  (k  X  E) 


}  . 


(35c) 


Since  j  =  a  •  E,  Eq.  (35)  may  be  rewritten  as: 

(o)^  -  c^k^)  I  •  E  +  c^k(k  •  E)  +  4::  i(jj  a  •  E  “  0  , 


(36) 
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which  may  be  expressed  as 


R  •  E  =  0  .  (37) 

R  is  called  the  dielectric  tensor  since  it  gives  the  dielectric  response  of  the 
plasma. 


In  Eqs.  (33)  through  (37)  E,  J,  and  v  are  the  first  order  in  perturbation  theory 
Fourier  coefficients  of  the  electric  field,  current  density,  and  velocity  of 
species  1  of  the  previous  section;  v^jj  represents  the  electrostatic  and  v(2)  the 
magnetic  contribution  to  the  perturbed  velocity.  For  simplicity  of  notation  we 
have  dropped  the  tilde  and  superscript  (1)  and  below  we  designate  a  generic  plasma 
species  by  the  subscript  I  Instead  of  1  as  In  the  previous  subsection.  Here  sub¬ 
script  i  Indicates  background  ions  unless  it  appears  as  a  summation  variable. 

To  solve  Eq.  (33)  we  evaluate  Eq.  (35),  substitute  the  result  into  Eq.  (34)  for 
each  plasma  species,  sum  over  species  and  finally  substitute  that  result  into 
Eq.  (33),  thereby  obtaining  an  equation  of  the  form  R  •  E  =  0  which  is  then  solved 
for  the  normal  frequencies  under  suitable  approximations. 


For  background  ions  and  electrons  the  streaming  velocity  is  zero  =  0,  where  I 
denotes  either  ions  or  electrons).  In  this  case  Eq.  (35c)  is  identically  zero. 
After  some  algebra,  the  right  hand  side  of  Eq.  (33)  for  each  background  species  is 
found  to  be 


-Awiuijj^ 


(38) 


4nN  e^  e  B 

In  Eq.  (38),  H  ( — jjj - )  is  the  plasma  frequency,  is  the 

1*  1 
cyclotron  frequency  and  we  have  introduced  Cartesian  coordinates,  A  similar 

calculation  for  the  streaming  ions,  neglecting  the  magnetic  terms  (we  perform 

simulations  using  an  electrostatic  code)  gives: 
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Here  the  subscript  b  refers  to  beam  ions.  Subscripts,  1,  e,  whenever  they 
appear,  designate  background  Ions  and  electrons,  respectively.  The  magnetic 
and  the  full  electromagnetic  contribution  to  the  current  density  have  also 
been  evaluated.  The  results  will  not  be  written  down  since  we  shall  subse¬ 
quently  be  Interested  In  the  electrostatic  approximation.  Dropping  the  mag¬ 
netic  terms  in  the  current  density  amounts  to  taking  kV/u)  «  1  as  a  comparison 
of  the  electrostatic  and  magnetic  terms  shows.  Making  use  of  Eqs.  (33),  (34), 
(38),  and  (39),  the  matrix  elements  of  the  dielectric  tensor,  as  defined  by 
Eq.  (37),  may  be  written  down  explicitly: 
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We  next  wish  to  solve  for  the  normal  frequencies.  To  obtain  tractable 
expressions,  we  look  for  normal  frequencies  w  such  that 


^i,b  ^e 


Making  use  of  Eq.  (41),  the  dielectric  tensor  becomes: 
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For  the  normal  frequencies  we  demand  det  R  »  0,  i.e., 

1^1  “  ^Scx^yy^zz  ”  ^xy®-yx^zz  ~  ^xzRzx^yy  =0  •  (43) 

To  obtain  a  reasonably  simple  dispersion  relation  we  further  restrict  our¬ 
selves  to  frequencies  co  such  that  co/kc  <<  1  (which  is  also  a  consequence  of 
the  electrostatic  approximation).  We  are  interested  in  beam  ions  nearly  per¬ 
pendicular  to  the  magnetic  field  as  this  produces  the  most  electron  heating 
and  so  we  take  kzAx  <<  1*  Carrying  out  the  algebra  in  Eq.  (43),  dividing  by 
c'^kx^k^  and  keeping  the  largest  contributions  only,  we  arrive  at 
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Noting  that  k,Av  <<  1  and  that  w  »  u)  ,  ,  Eq.  (36)  may  be  given  approximately 
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This  result  is  to  be  compared  to  the  result  of  McBride  and  coworkers  (1972). 
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The  difference  between  Eqs.  (45)  and  (46)  is  (k.2/kx)^  versus  (kg/k)^,  which  is 
negligible  because  k^  "  k^^.  Whereas  the  McBride  result  is  obtained  as  the 
cold  limit  of  a  warm  theory,  the  derivation  here  has  as  starting  point  a  cold 
plasma.  Note  also  that  in  McBride's  work  there  are  no  background  ions. 


Let  us  assume  equal  masses  and  charges  for  streaming  and  background  ions,  and 
equal  charges  for  electrons  and  ions.  Furthermore,  let  us  define  the  lower 
hybrid  frequency  by 
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Upon  dividing  Eq.  (45)  through  by  1  +  Upe^/fle^»  using  Eq.  (47) 

=  (m, /m  )(N  /N,  )ai^,  and  oi^,  =  (N. /N,  )(d^  ,  Eq.  (45)  may  be 

pe  1  e  e  1  pi  pb  b  1  pi  ’  ^  ■' 


and  noting  that 
rewritten  as 


(1  +  6^)  ^  +  Sj. 


(oj  -  k  V)^ 
X 


(48a) 


where  6^  =  (k2/kjj)^(mi/rag)  (Ng/Nj^)  (48b) 

Gr  =  N^/Nj  ,  (48c) 

This  result  is  to  be  compared  to  Abe  and  Machida  (1985).  The  difference  is 
that  there 


(49) 


For  kg  «  kjj,  Eqs.  (47)  and  (49)  are  essentially  identical. 


The  dispersion  relation  (Eq.  (48))  is  easily  written  as  a  quartlc  polynomial  in 
a)(k)  and  solved  numerically.  A  plot  of  oj  versus  k^V  appears  in  Fig.  7.  Our 
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choice  of  parameter  normalizations  are  such  that  the  velocity  of  the  beam,  the 
mass,  cyclotron  frequency,  and  the  charge  of  the  electron  are  all  unity.  In 
addition  the  number  density  of  electrons  and  ions  Is  taken  to  be  the  same. 

The  figure  shows  the  real  (ai^)  and  Imaginary  (y)  parts  of  a)(k)  as  a  function 
of  kxV. 

Physically,  y  gives  the  unstable  growth  rate  of  electrostatic  waves.  We  are 
Interested  in  domains  of  k  space  where  y  is  non-zero  for  the  following  reason. 
In  order  that  the  electrons  be  heated  as  required  by  the  CIV  hypothesis,  the 
plasma  must  be  able  to  transfer  energy  from  the  beam  to  the  background 
electrons.  This  Is  achieved  via  a  plasma  Instability  which  pumps  energy  from 
the  beam  to  the  electrostatic  field  which  is  then  Imparted  to  the  electrons. 

In  Fig.  7  the  maximum  growth  rate  Is  y  =  0.02  and  is  achieved  at  kjjV  =  0.095. 
From  the  well  known  relation  from  finite  Fourier  transforms  k  =  (2n/L)n,  with 
n  an  Integer,  we  can  readily  calculate  the  period  of  the  physical  length  of 
the  system  which  leads  to  Instability, 

3.3  COMPUTER  SIMULATION  OF  ION  BEAM  STREAMING  INTO  A  BACKGROUND 
ION-ELECTRON  PLASMA 

During  the  year  we  modified  the  ESI  plasma  simulation  code.  This  code  was 
developed  by  Blrdsall  and  coworkers  (Blrdsall  and  Langdon,  1985)  for  partlcle- 
In-cell  simulations  of  plasma  effects.  The  code  Is  one-dimensional  in  posi¬ 
tion  space  and  was  originally  two-dimensional  in  velocity  space  with  the  mag¬ 
netic  field  directed  perpendicularly  to  the  plane  of  the  velocity  vectors.  We 
have  modified  ESI  to  be  three-dimensional  in  velocity  space,  so  that  we  now 
consider  four  of  the  six  dimensions  of  phase  space  for  each  particle.  We  now 
allow  the  magnetic  field  to  be  at  any  angle  with  the  ion  beam  velocity  vector; 
we  take  the  ion  beam  velocity  vector  to  be  along  the  x  axis,  where  x  is  the 
one  position  variable  in  the  code.  With  the  modification  we  can  now  simulate 
the  conditions  that  we  developed  the  linear  theory  for  in  Subsection  3.2 
above.  We  ignore  collisions  with  neutral  molecules  for  the  results  in  this 
section.  Such  collisions  are  considered  in  Subsection  3.5. 
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The  case  of  the  fast  ton  beam  streaming  into  a  plasma  represents  a  version  of 
the  CIV  experiments  where  a  fast  neutral  beam  is  injected  into  the  ambient 
plasma  of  the  ionosphere  about  the  space  shuttle.  This  is  because  the 
ionization  of  the  fast  neutrals  produces  a  beam  of  fast  ions,  as  there  Is 
little  momentum  change  of  the  neutral  molecule  when  it  is  ionized.  For  the 
simulation  we  take  the  electron  plasma  frequency  a>p  of  the  ambient  plasma  to 
be  equal  to  the  electron  cyclotron  frequency,  f2g.  For  example,  if  the 
magnetic  field  B  is  0.33  Gauss  then  the  ambient  electron  density  n^  is 
1.06  X  10^  cm3,  or  n^  is  1.98  x  10^  cm3  for  a  B  of  0.45  Gauss,  as 

(ijp/ng  =  3.21  X  10-3  nel/2/B 
when  B  is  in  Gauss  and  n^  is  in  number/cm3. 

McBride  and  coworkers  (1972)  show  that  the  maximum  growth  rate  y  occurs  when 
the  angle  iji  between  the  beam  velocity  vector  and  the  magnetic  field  is 

'^'max  “  90“  -  sin-^(rae/mi)^/2  ^ 

where  m^  and  mj  use  the  masses  of  the  electron  and  the  ion,  respectively.  We 
now  present  the  results  of  simulations  for  two  values  of  ip.  We  use  an  arti¬ 
ficially  high  mg/ra£  ratio  of  0.01,  which  speeds  the  computation  while  still 
preserving  the  essential  physics  that  ions  are  much  more  massive  than  elec¬ 
trons.  Also,  we  take  the  ions  to  be  unmagnetized  (i.e.,  they  have  an  infinite 
gyroradius),  which  correctly  represents  the  other  important  physical  differ¬ 
ence  between  the  ions  and  electrons.  The  maximum  growth  then  occurs  at 
'('max  “  84.3“.  We  discuss  results  for  ip  values  of  3  and  87  deg  -  one  value 
slightly  greater  than  ipmax  value  that  is  nearly  perpendicular  to  ipmax* 

In  both  simulations  oip  =  Qq,  the  velocity  of  the  fast  ion  beam  is  100  times 
the  mean  thermal  velocity  of  the  ambient  Ions  (which  Is  equal  to  the  mean 
thermal  velocity  of  the  electrons  so  that  the  electron  temperature  is  less 
than  the  Ion  temperature,  Tg  >»  0.01  T^),  and  the  number  density  of  the  fast 
ion  beam  is  10  percent  of  that  of  the  ambient  ions. 
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Figure  8  shows  the  growth  of  the  electric  field  energy  for  the  two  angles. 

This  shows  that  the  field  energy  builds  up  to  a  value  that  is  several  times 
larger  for  ij;  =  87  deg  than  it  does  for  ip  =  3®.  We  can  compare  the  logarithmic 
growth  rate  with  that  predicted  by  y  in  Fig.  7  for  =  87°.  The  predicted 
value  from  linear  theory  has  y  =  0.02  Oq  while  the  slope  of  Fig.  8a  (which  is 
related  to  y  by  the  factor  of  2.303/2)  gives  y  =  0.018  Og  for  times  for  times 
from  110  to  l^O/fig.  At  later  times  from  140  to  180/0e  the  slope  is  larger  as 
non-linear  effects  dominate;  here  y  =  0.038  Og  for  tp  =  3®,  y  =  0.012  Og  for 
times  from  115  to  140/0g  and  y  =  0.034  ftg  for  times  from  150  to  180/0g.  Thus 
the  field  energy  builds  up  even  for  ^  far  away  from  tj/max’ 

When  the  field  energy  increases  the  ion  beam  develops  wavelike  structures  as  a 
result  of  the  electrostatic  field.  The  field  then  begins  to  trap  particles, 
and  the  growth  saturates.  Figure  9  shows  the  velocity  along  the  direction  of 
the  fast  ion  beam  for  i|/  ®  3®.  The  velocity  of  each  particle  is  plotted  as  a 
point,  so  the  coherent  effects  are  evident.  The  wave  structure  at  v^  “  1  is 
in  the  fast  ion  beam,  while  the  ambient  ions  have  nearly  zero  velocity.  The 
electrons  also  have  a  wavelike  structure  in  their  velocities. 

Figure  10  shows  a  comparison  of  the  electron  kinetic  energies  for  the  two 
angles.  Here  the  perpendicular  motion  gives  over  twice  as  much  energy  to  the 
electrons,  but  it  is  interesting  that  the  nearly  parallel  motion  can  still 
heat  the  electrons  as  much  as  it  does.  The  Increase  for  times  greater  than 
350/0g  is  caused  by  a  plasma  effect  that  transfers  energy  back  and  forth 
between  the  beam  ions  and  the  electrons.  Additional  calculations  show  that 
the  kinetic  energy  of  the  electrons  has  a  maximum  of  approximately  102  at  a 
time  of  450/0g  and  it  then  declines  to  a  value  of  68  at  a  time  of  530/0g  (the 
energy  units  are  normalized  so  that  the  initial  energy  in  the  ion  beam  is 
approximately  349). 

The  increase  in  kinetic  energy  of  the  electrons  occurs  at  the  same  time 
as  the  drift  energy  of  the  fast  ions  decreases,  which  is  shown  in  Fig.  11. 

Note  that  the  initial  drift  energy  is  the  same  for  both  angles.  Only  a 
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Figure  8a.  Growth  of  the  electric  field  energy  at  two  angles  between  the 
magnetic  field  and  the  beam  velocity.  Time  Is  measured  In  units 
of  the  reciprocal  of  the  electron  cyclotron  frequency  Dg,  and  the 
energy  units  are  such  that  the  Initial  drift  energy  of  the  fast 
Ion  beam  has  a  value  of  349.  Figure  8a  Is  for  -  87“ 
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Figure  8b.  Growth  of  the  electric  flel-^  ..gy  at  two  angles  between  the 
magnetic  field  and  the  beam  velocity.  Time  Is  measured  In  units 
of  the  reciprocal  of  the  electron  cyclotron  frequency  and  the 
energy  units  are  such  that  the  initial  drift  energy  of  the  fast 
ion  beam  has  a  value  of  349.  Figure  8b  Is  for  3® 
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Figure  9.  Velocity  along  the  direction  of  the  fast  ion  beam  as  a  function  of 
distance  at  a  time  of  199.975/^0.  The  fast  ions  initially  all  had 
Vx  =  1.0;  now  their  values  have  a  coherent  structure.  The 
electrons  initially  had  a  distribution  of  velocity  with  a  mean 
value  of  0.01  -  similar  to  that  of  the  ambient  ions 
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Figure  10b.  Electron  kinetic  energy  at  two  angles  ip  between  B  and  Vx.  T; 

time  units  are  l/iig;  energy  can  be  normalized  to  the  initial 
drift  energy  of  the  ion  beam  by  dividing  by  349.  Figure  10b 
is  for  ■  3“ 
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Figure  11a.  Ion  drift  energy  for  two  angles  between  B  and  v^.  The  time  t 
Is  measured  In  units  of  the  reciprocal  electron  cyclotron 
frequency  Figure  11a  Is  for  <1/  =  87® 
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Figure  11b.  Ion  drift  energy  for  two  angles  tp  between  B  and  Vx*  The  time  t 
Is  measured  In  units  of  the  reciprocal  electron  cyclotron 
frequency  £1-.  Figure  11b  Is  for  ili  ■  3". 


fraction  of  this  energy  is  removed  from  the  beam.  This  fraction  is  smaller 
for  =  3“.  Additional  calculations  show  that  the  drift  energy  decreases  to 
112  and  then  increases  to  215  and  decreases  to  130  (where  349  is  the  initial 
beam  energy)  at  times  between  400  and  800/Qg.  In  these  simulations,  the  total 
energy  is  conserved  to  about  0.1  percent. 

These  results  illustrate  some  of  the  detail  available  from  the  simulations. 
Note  that  there  is  electron  heating  over  a  wide  range  of  angles.  However, 
the  electron  energy  gain  is  smaller  for  <|/  =  3®  than  for  =  87®.  Note  also 
that  there  is  a  delay  of  some  170/Qg  before  the  electrons  heat  appreciably. 
Then  the  electrons  are  quickly  heated  in  a  time  of  about  50/Qg. 

3.4  WAVES  IN  A  "HOT"  PLASMA  IN  A  MAGNETIC  FIELD 

Following  the  philosophy  for  waves  in  a  cold  plasma,  we  study  the  properties 
of  small  amplitude  disturbances  propagated  through  a  "hot"  plasma  in  a  uniform 
magnetic  field.  A  uniform  unbounded  plasma  is  assumed.  Our  approach  is  to 
solve  perturbatively  the  relativistic  Vlasov  kinetic  (collisionless)  equation 
assuming  probability  conservation  and  the  relativistic  Maxwell's  equations. 

We  begin  with  a  brief  review  of  the  general  problem  in  order  to  establish  the 
notation  and  the  geometry  and  then  specialize  to  the  system  of  interest  for 
the  CIV  experiment. 

3.4.1  General  Theory 

The  relativistic  equation  of  motion  for  a  particle  of  charge  e  and  rest  mass 
mo  is 


dp 

Jt 


where 


=  e(E  +  i  V  X  B)  s  G 
Ttn^  V  ;  Y  =  [l  -  (v/c)^l 


-1/2 


(50a) 

(50b) 
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If  w<3  let  the  probable  number  of  particles  of  a  given  species  in  the  phase 
space  volume  element  dx  dp  at  t  be  given  by 

dN  =  n^f(x,p,t)  dx  dp  (52) 

then  the  distribution  function  f  satisfies  a  relativistic  generalization  of 
the  collisionless  Vlasov  equation,  namely 


In  o;:der  to  solve  the  complex  coupled  system  of  differential  Eqs.  (53) 
and  '54)  ve  attempt  a  perturbative  solution  about  a  homogeneous  plasma  with  a 
uniform  magnetic  field  Sq,  no  electric  fields  present,  the  zero  order  dis¬ 
tributions  have  no  angular  dependence,  and  the  total  current  carried  by  the 
distributions  must  vanish  in  the  zero  order  system.  Letting 
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(56) 


f  =  fO  . 

E  = 

B  =  B  +  B^^^ 
o 


and  substituting  f,  E,  B  into  the  system  (53)  and  (54),  we  obtain  a  system  of 
linear  equations  which  can  be  solved  using  Fourier-Laplace  transforms,  defined 
by: 


(E,B,J,f)  =  J  dt  e"®^  j 


dx  -iR*x  r=(l)  5(1)  x(l)  ,(1) 


(E<^>,B(l>,j(l>,f(l>) 


with  R  (s)  >  0,  and  the  inversion  formula 

gV  / 

+®  o+i“ 

=  J  dfc  J  ^  e®^E,B,J,f)  .  (58) 

_co  a-i® 

From  here  on  we  deal  only  with  transforms,  so  there  should  be  no  confusion 
between  the  total  fields  in  Eq.  (56)  and  the  quantities  on  the  left  of 
Eq.  (57).  Also,  unlabeled  suras  mean  summation  over  charged  species.  The 
motivation  of  solution  is  the  same  as  the  cold  case,  that  is  eliminate  all 
variables  except  one  (the  electric  field  is  a  good  choice),  write  the 
resulting  equations  as  an  eigenvalue  problem  and  solve  for  the  eigenvalues  and 
eigenvectors.  The  result  of  that  analysis  is 


2  +® 


(s^  +  c\^)  E  -  c^(lc*E)E  +  Alls  ^  J  dp  p  J  d+'G(^',+) 


[.  , 

E  -  -  V'  X  (R  X  E)  •  — -  =  I  ,  (59a) 

®  J  3p' 


where 
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I»si  +  icRxB  +  Alls  £  ^  ^  J  dp  p  J  d+'G(+',+) 

o  o 


+“ 


S 


3f 

-  —  v'  X  B  •  — ^ 
3p'. 


In  Eqs.  (59a)  and  (59b) 


;{ 


Q 


G(^,*')  =  e 
with  properties 

r  G(+-»*'  -»-»)->  0  , 

1  G(+-»+',^)  -»  1 


k|V| 

(♦-♦') - CT  (sin  ♦'  -  sin 


♦)| 


All  finite  ♦ 


e|B  I  a. 

Q  ■  - °  ■  -2  is  the  relativistic  cyclotron  frequency, 

"o^  Y 


and 


+  ® 


(i,B,i)  ■  J  [E^^^t=0),B^^>(t=0),f^^^t=0)) 


(2ll)' 

B. 


=  (0,0, Bq) 

k  =  (k^,0,k^) 


A-6696 

Figure  12.  Coordinate  system  with  magnetic  field  along  the  z-axis. 


(59b) 


(60a) 


(60b) 


(61) 


AO 


Note  that  the  initial  values  e,B,g  have  been  collected  in  I.  In  addition  we 
have  adopted  a  coordinate  system  with  the  constant  magnetic  field  along  the 
z-axis  and  It  in  the  x-z  plane. 

Equation  (59)  is  conveniently  written  in  the  form 

(s^  +  c^k.^)E  -  c^It(It*E)  +  4ns  o*E  =  I  (62) 

where  o  is  the  conductivity  tensor,  defined  through 

J  =  (63) 

Equation  (62)  is  also  in  the  form 

J-E  =  I  (64) 

where  ^  is  a  dyadic  which  may  be  written  as  a  3x3  matrix,  if  E  and  I  are 
written  as  three-component  column  vectors.  Equation  (64)  has  solution 


where  R  is  the  determinant  of  the  matrix  and  R^^  is  the  transpose  of  a 
matrix  whose  elements  are  the  co-factors  of  the  corresponding  elements  of  Rjj. 
The  elements  of  Rjj  are  to  be  obtained  from  Eqs.  (62)  and  (63).  The  last  term 
on  the  Ihs  of  Eq.  (62),  i.e.,  Eq.  (63)  is  the  complicated  part. 

Subsequently  we  will  be  interested  in  implementing  the  electrostatic 
approximation  7  x  E  »  0,  which  implies  the  relation  E  x  E  =  0  for  the  Fourier 
coefficients.  Assuming  a  zero  order  distribution  fo  =  fo(P|,P||),  Eq.  (64) 
can  be  evaluated  explicitly.  The  result  is  ~ 


t 

R 

R 

XX 

xy 

xz 

R 

R 

Vx 

yy 

yz 

R 

R 

zx 

zy 

zz 

where 
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„  2  2 
R  =  s  +  c  i 

XX 


h  -■’III 


n^aVoyapj) 


ik, ,V| I  +  inQ) 


nQ^J  J'(0f^/8p,) 
n  n'  0  *^1 


R  =  -R 
xy  yx 


■  I  ^  -"ll  J  ■’i  -”1  *  lk||V||  »  iHS 


R  =  -c  k,  k,  ,  -  2 IIS 
xz 


\  "d  \  r  ,  r  2  .  neVofQ/3p||) 


O  2  2 

R  =  s  +  c 

yy 


2  ®  ®  2(j;)^Ofo/3p,) 

. •T-.i 


-00  0 


V*  V  r  p  o  J  ^  ®P  I  I  ) 

I  I  J  -■’ll  I  "i  -"i  is  /ik,,v---rra 


n2V(3fQ/3P|) 


R  =  -c 
zx 


'ki|ki  2iis  ^  qP  ^  J  P||dPj|  J  P|  k|V|(s  +  ik|  |V|  j  +  InB 

n  ~  ^ 

^  -  -  v;<v>p|) 

"sy  ■  2  ^  2  i  "I  I-"!  I  i  "l  -"I  (S  *  *  in5) 


n  2  2 

R  =  s  +  c 
zz 


2  V  “d  Y  r  .  r  . 

'‘|-'’“2s^  zJ  ‘’ii-'’ii  J  "1  -"1 


The  above  sums  are  defined  for  R^s  >  0,  and  the  argument  of  the  Bessel 
functions  is  k|V|/a  =  k|P|/m^a^.  The  unlabeled  sums  run  over  charge  species. 
In  arriving  at“Eq.  (67)7  use  has  been  made  of  the  Bessel  function  identities 
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sin 


2n 

J  ^-lz[sin(«fr-a)  -  sin  +1 

0 


sin 

cos 

cos 


<|>  sin(<^-a) 
♦  cos(+-a) 
Ip  sin(4>-a) 
<j>  cos  (♦-a) 
1 

sin  <|> 
cos  <J> 

sin(<j>-a) 

cos(  <J>-a) 


+  CO 


'"I 


ina 


n=-® 


(j;)- 


j  j' 
n  n 

J  J' 
n  n 


in 
z~ 
in 

n  ,2 
~2  n 

-iJ  J' 
n_n 
n  j2 

z  n 

iJ  J' 
n  n 

n  j2 

z  n 


(68) 


The  equation  R=0  will  be  studied  below  and  will  be  described  as  the  t->® 
"dispersion  relation",  but  we  must  bear  in  mind  that  the  residues  from  the 
poles  of  the  reciprocal  of  R  give  only  a  part  of  the  asymptotic  t-^  solution. 
These  results  reduce  to  the  "cold"  wave  results  in  the  double  limit  v^l  and 
T-K).  The  first  limit  is  simply  the  nonrelativistic  limit.  For  a  more 
detailed  derivation  of  these  results  the  reader  should  consult  Montgomery  and 
Tidman  (1964). 


3.4.2  Dispersion  Relation  for  a  "Hot"  Plasma  in  a  Magnetic  Field 

The  system  consists  of  a  background  electron-ion  plasma  and  streaming  ions. 

We  fix  coordinates  on  the  background  ion-electron  plasma  and  take 

B  =  (0,0,Bq)  =  BqZ 

R  =  (k^,0,k^)  =  k|X  +  k||Z  (69) 

u  =  (V,0,0)  =  V|X 

In  this  calculation  we  make  the  electrostatic  approximation.  This  implies 
that  only  the  terms  (67)  are  present  and  «/kc  <<  1.  More  explicitly  we  look 
at  frequencies 
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kV,  «  «  <<  kc  , 


(70) 


where  Vj  refers  to  the  streaming  velocity.  If  the  initial  perturbations  e,B,g 
are  identically  zero,  then  I  s  0. 

The  background  distribution  is  taken  to  be: 


^VkgT,,  2m^Yk3T, 


2nmoY  kgT|  [2nm^Y  kgT|  | 


(71a) 


while  the  streaming  ion  distribution  is: 


T||^ 

T,  -K) 


->  5(P| |)6(P|-P) 


(71b) 


The  delta  function  normalization  is 


I  dPiiMPii)  - 1 


2ii  ® 


J  d*  J  dpj  P|  5(P|-P)  =  1 


0  0 

Letting  |R|  =  0  and  keeping  only  the  largest  contributions  (recall  kc  >>  w)  in 

4  2  2  2 

analogy  to  the  "cold"  plasma  calculation  we  find,  upon  dividing  by  c  k.k  s  : 


/  if  /  u  \ 

k  .  til  R  .  ^  1  [b 

XX  k,  zz  k,  s  ( 


+  R 

xz  zx 


)  =0 


where 
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"xK  ■  1  2  r  2  J  ''"ll  J  ‘*T 


2  3  2 

n^Q-'j^oyap.) 


-  -  k|  V|(s  +  ik,,,v,,  +  inS) 


;  -  1  2n  V  “p  X  r  .  r  . 

*^zz  "  s  z  2  J  ‘*^11^11  J  (s  +  1^1  I  + 


nsVof  zap, ,) 


2n)/)  [  dp,,  f  dp.pf 

L  %  L  J  II  J  i  1  k|V.(s  +  ik. ,v. .  +  in 


/  +®  ® 

H  I  I  '*'’ii'’ii  1  Vi 


ns^j^(af  zap.) 

n'  o 

I V  I (s  +  ik| I V  I  I  +  inQ) 


We  further  restrict  our  attention  to  the  case  k||Zkj^  «  1  and  keep  lowest 
order  terms  only. 


[s  +  inQ  +  I'^l 


1  ‘Nrii, 

1  —  — — . 1....-^  .f  * 

s  +  inQ 


Under  these  circumstances,  the  matrix  elements  R  in  Eq.  (74)  simplify. 

Since  fo  is  even  in  pj j  for  either  of  the  distributions  (71a, 71b),  it  is  clear 


that  R  =  R  =  0. 
xz  zx 


The  dispersion  relation  becomes; 


R  + 

XX 


(k 

•I  R 
U  2 

ri 


The  matrix  elements  R^^,  R^^  can  be  found  by  evaluating  the  momentum  integrals 
and  simplifying  the  sums  over  n.  We  get,  for  background  species 
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00 


R  =1  +  2 

XX 


2  n2 - ’T- 

w  m  Q 

poo  m  Q 

-JT - 1  ®  °  ° 

rVl'^i 


2t 

n  I 


n  i 

5 


V 


n=l 


R  =  1  + 
zz  s 


2  m  Q 

we  o  o 
P 


2s 


5?^ 


n=l 


f  2  2^ 

v^k  T  K 
Y  Kgl.K, 


where  the  argument  of  I  is  I  =  I 
**  n  n  n 


T 

m  S 
o  o 


(77a) 


(77b) 


In  obtaining  Eq.  (77),  we  have  noted  that 


(i) 


+  ® 


J  ■‘■’I 


2"„r  Vi 


2niii^Y  kjT. 


(78a) 


(li)  J  dP|P|j2 


^"o''  V|  Z. 

'  -  -  "o''  ''b^I  ' 


7  ?T 

m  Q 
\  o  oj 


(78b) 


see  Watson  (1966,  pg.  395) 


with 


f  2  2'\ 

rVlV 


I  =  I 
n  n 


— r 

m  Q 
o  o 
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(iii) 


Jt 

r^nl  rinl 

5  s  +  'in-a  =  ®  ^  ^  2  T ^ 

^  ^s+nQ  ^s+n 


r-i  n  I 

2  ^ ^ 

^  s  +  n 


The  streaming  ions  are  taken  to  be  cold  and  unmagnetized,  i.e., 

B  =  0  =>  Vj/Q  -»  ®,  i.e.,  an  infinite  Larmor  radius.  Rather  than  proceed  as 
above  we  present  a  derivation  of  the  weak  magnetic  field  limit  of  Eq.  (59), 
and  then  evaluate  Rxx»  ^zz  distribution  of  Eq.  (71b).  In  the  weak 

magnetic  field  limit  we  can  expand  the  third  term  of  Eq.  (59)  by  successive 
integrations  by  parts.  Thus  suppose  we  write 


E  -  i  V'  X  (R  X  E)  •  — 

®  8p' 

The  last  term  on  the  left  of  Eq.  (59)  is  then 


5  -n 


dp  -  I 
Y 


(81a) 


where 


r  .r’mnr).*  f 

>*  ^  (s  +  ik-v') 


A(<|)')  ac 


(81b) 


In  arriving  at  this  last  result  we  have  used  Eq.  (60)  to  obtain 


^  (s  +  ik* v' ) 


A7 


If  we  now  successively  integrate  Eq.  (81b)  by  parts  and  substitute  G/Q  using 
Eq.  (60),  a  series  results,  namely: 


I  =  + 


A(*) 


s  +  ilc-v  (s  +  ilc'v) 


3 


(  \ 

A(+) 


U  +  ilc’v  j 


(s  +  ilc'v)  ^ 


(s  +  ilc'v) 


A(*) 


ts  +  ilc'V 


(83) 


where  we  also  have  noted  G( ♦=♦';♦')  =  1* 


To  zero  order  in  Q,  Eq.  (59)  reduces  to 


/  2  2.2,=  2 

(s  +  c  k  )E  -  c 


Ic(Ic’E)  -  s  ^  f  - -  [e  -  -  V  X  (It  X  E)]  •  — ^  =  I 

^  P  **  y(s  +  ilc'v)  •-  s  ■*  3p 


For  kv  «  0),  Eq.  (84)  becomes: 


(s^  +  c^k^)E  -  c^Ic(lc’E)  -  s 


5  r  df  — E — _  E  •  . 

^  ^  y(s  +  iE*v)  3p 


y(s  +  ik*v) 

Thus,  using  cylindrical  coordinates,  we  thus  must  evaluate: 


I  (85) 


^2  +”  ®  2ii 

'‘P||  J  J  '‘♦P| 


2a 
cos  9 


(3f^/3P|) 
s  +  ik.v, 


2  +® 


2ii 


V’Wp  P  P  (3f/9Pl|) 

“xz  ■  2  I  -‘Pi  I  1  “PiPi  1  o+Pi  I  X !  ik|;| 


(86a) 


(86b) 


0 


0 


where  is  given  by  Eq.  (81b)  subject  to  the  normalization  Eq.  (82). 
Straightforward  integrations  lead  to 
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(87a) 


R  =  + 

XX 


^  (1)  2 
)  s 

^  ^  (s  +  ik.v,)' 


R  =  + 
zz 


L 

\  _P  s 
^  ^  (s  +  ik.v.) 


(87b) 


where  V|  =  P|/Yino.  Subsequently,  we  shall  be  interested  in  the  nonrelativistic 
limit  Y=l'  ~ 

Restricting  ourselves  to  frequencies 


^i j b  ^  ®e 


and  writing 


s  =  -iw 


we  obtain  for  the  physical  system  under  consideration 


o  2 

2w  m 

R  =  1  .  -  -P-g-  % 

k  T  k^ 

Vie*"! 
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Vli"! 


2 

(li  .  e 
pi 


m.QT  +I  .  „ 

11  01  2 

~T~  ^  1 
<«>  u 


Wt  0)  -  K  ,  V  I 


(90b) 


In  Eqs.  (90a)  and  (90b)  the  subscripts  e,i,b  refer  to  electrons,  background 
ions  and  beam  ions,  respectively. 

For  the  sake  of  compactness  of  notation  we  introduce  a  number  of  definitions: 


2  -  _ ^  ®  f  \ 

2m  Q  -2  ^  „ 

,  o  m  Q  \  ^  2 

^  kTIc  A  ^  ' 

2  B  ,  ll/n 

bJ 


te  m 


(electron  thermal  velocity) 


(92a) 


rtf  T 
m. 

V  1  2 


(ion  thermal  velocity) 


(92b) 


a  - ^ -  (lower  hybrid  frequency) 

« 

I  +  P®  Z 
e 

9  «  m.  N 

and  6  a  (k.  .  /k.  )  —  =5^ 

II  i  "‘e  ”i 

In  addition  recall  the  following  plasma  frequency  relations: 


2  "’i  ^e  2 

*^pe  ~  m  N.  *^pi 

2  "’i  ^b  2  2 

w.=  —  ae«. 

pb  mj^  pi  r  pi 


(95a) 


(95b) 
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Substituting  Eqs.  (90a)  and  (90b)  in  the  dispersion  relation  (Eq.  (86)), 
dropping  (k||/k|)^  as  well  as  terms  of  order  (k| |/kj we  obtain  upon 
making  use  of  difinitions  (91)  through  (95)  our  finil  simplified  dispersion 
relation 


(I) 


Ih 


*  s  e 


'  I 


oe 
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Introducing  the  shorthand  notation 


♦  =  ». 
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Z,i  *  e  ^ 


e  ) 


■•Oe 


(96) 


(97) 


6  a  1  +  8' 


2  “ih 

T 

e 


^3e 


The  dispersion  relation  (Eq.  (97)  may  be  rewritten  more  simply  as; 


(98) 


c  w,. 

r  Ih 


~2  * 

M  (w  -  k,V|  )' 


=  A 


(99) 


For  computational  purposes  Eq.  (99)  can  be  expressed  as  a  quartic  polynomial: 


iJ*  -  2Ak|V|»^  +  [a(k|V|)^  -  ♦  -  e^wyo)^  +  2i  kj^v^w  -  Kk^v^)^  =  0 


(100) 
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3.4.3  Parametric  Study  of  CIV 


In  this  section  we  study  the  "warm"  plasma  dispersion  relation  developed 

previously.  From  an  inspection  of  the  dispersion  relation,  the  independent 

parameters  are  k.|,/k.,,  m  /m.,  e  ,  «  /Q  ,  v^,  /v^.  . ,  and  v,/v^,  In  what 
I  I  I  e  1  r’  pe  e  the  thi’  |  thi 

follows  we  examine  tHe  dependence  of  the  maximum  growth  of^uns table  waves  on 
each  of  these  parameters. 

The  electron  densities  and  magnetic  field  lie  in  the  ranges 
10^  <  ne  <  5  X  10^  #/cm3 

0.3  <  B  <  0.45  gauss 


so  that 


0.7  <  Wpe/Qe  <  24 

Adopting  the  nominal  values  of  n  =  5  x  10^  cm“^  and  B  =  0.38  gauss  gives  us  a 
nominal  (Ope/^Jg  =  6.  The  gas  species  that  are  to  be  ejected  in  the  experiment 
are  Ne,  Xe,  NO,  and  CO2.  Assuming  the  plasma  ions  are  essentially  all  O'*'  gives 
us  a  range  for  the  ratio  of  beam  ion  mass  to  background  plasma  ion  mass  of 

1  <  m. /m.  <  8 
-  b  1  ~ 

This  ratio  appears  in  Sj-  =  (Nb/Nj ) •  (mj/mj,) ,  which  determines  a  scaling  law, 
namely  a  given  tj-  may  be  obtained  by  increasing  N^/Ni  and  decreasing 
m^/m^^  or  vice-versa.  Taking  the  beam  velocity  V|=l  unit  gives  a  range  for  the 
ion  and  electron  thermal  velocities 

0.07  <  v^,  ,  <  0.09 
-  thi  ~ 

1.2  <  V,.  <2 

~  the  ~ 
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Our  nominal  values  are  v^jji  =  0.088,  and  =  1.4.  The  thermal  velocities 

were  calculated  from  v^^e  =  x  10^  gmj  =  9.79  x  10^  p-1/2 

where  Tg,  Tj  are  the  electron  and  ion  temperatures  in  eV  and  p  is  the 
ion  mass  expressed  in  units  of  the  proton  mass.  The  range  of  temperatures  of 
interest  are  0.08  <  Tg  <  0.2  eV,  0.05  <  Tj  <  0.10  eV.  Normalizing  thermal 
velocities  to  the  beam  velocity  and  using  a  model  mass  ratio  of  mg/m^  =  1/400, 
so  that  effectively  we  are  using  an  electron  mass  that  is  137  times  as  large 
as  the  physical  mass  of  the  electron,  gives  the  above  range.  We  summarize 
first  the  results  for  the  dependence  of  the  maximum  growth  rate  on  e,-.  The 
fixed  parameters  are 

k  /k  =  0.05  =  (m  /m.)^^^  V,=  1 

z  X  '  e  i'  I 


me/mi  =1/400 

“pe'^^e  ~  ^ 
mi/mij  =  1/3 


vthi  =  0.088 

Vfhe  =  1-^ 


The  results  are  summarized  in  Figure  13.  Ymax  increases  with  rapidly  at 
first  then  saturates  to  Ymax  =  0.018.  This  corresponds  to  Njj/Nj  =  3,  i.e., 
our  beam  density  is  three  times  as  large  as  the  background  density.  However 
the  growth  is  appreciable  even  for  low  values  of  Nij/Nj.  Figure  14  shows  the 
dependence  of  growth  rate  on  k||/k|  which  gives  the  inclination 
of  the  k  vector  relative  to  the  direction  of  particle  displacement  and  the 
self-consistent  electric  field,  i.e.,  the  x-direction.  In  the  range  of 
validity  of  the  dispersion  relation  i>  so  that  the  maximum  growth 

rate  increases  as  we  fall  away  from  the  perpendicular  direction.  However,  from 
preliminary  simulations  electron  heating  is  maximized  at  k|j/kj  = 
m^/m,^'^^.  Finally,  Figure  15  illustrates  the  dependence  of  the  maximum 
growth  rate  on  cOpg/Qg.  Here  again  it  increases  and  saturates  at  approximately 
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Figure  14. 


Parametric  dependence  of  maximum  growth  rate  on  ratio  of  beam 
to  background  number  densities. 


Parametric  dependence  of  maximum  growth  rate  on  ratio  of  parallel 
to  perpendicular  components  of  wave  vector. 


54 


Figure  15.  Parametric  dependence  of  maximum  growth  rate  on  ratio  of  electron 
plasma  to  cyclotron  frequencies. 

3.5  COMPUTER  SIMULATION  OF  MOLECULAR  AND  PLASMA  EFFECTS 

3.5.1  Introduction 


There  is  a  rich  variety  of  molecular  physics  involved  in  a  plasma  that  is 
producing  ionization  via  the  CIV  effect.  The  CIV  effect  seems  to  be  more 
complicated  than  first  thought.  For  example,  the  critical  velocity  for  a 
discharge  to  occur  is  larger  than  Vq  because  only  a  fraction  of  the  kinetic 
energy  of  the  fast  ions  is  transferred  to  the  electrons  (McBride,  et  al., 
1972;  Galeev,  1981)  and  because  some  of  the  electron  energy  goes  to  produce 
excitation  rather  than  ionization  (Newell  and  Torbert,  1985).  The  latter 
result  could  have  been  anticipated  from  the  results  of  radiation  physics 
studies  (e.g.,  Platzman,  1961).  However,  some  of  the  energy  that  goes  into 
excitation  may  be  recovered  as  a  result  of  ionization  of  metastable  states 
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(McNeil,  Lai,  and  Murad,  1987).  Another  factor  that  can  influence  ionization 
by  CIV  is  the  loss  of  ions  through  recombination  (Murad,  et  al.,  1986)  and  ion 
diffusion  (Caledonia,  et  al.,  1987). 


There  are  many  effects  that  must  be  considered  in  order  to  predict  CIV 
effects.  These  effects  include: 

1.  Seed  ionization  to  initially  produce  the  fast  ions  from  ionization  of 
fast  neutrals  in  an  event  with  low  momentum  transfer.  This  ionization 
could  come  from  photoionization,  ionization  by  fast  electrons  present  in 
the  ambient  ionosphere,  charge  exchange  collisions  with  ambient  ions,  and 
collisions  with  fast  neutral  molecules  or  atoms  from  the  ambient.  A 
comparison  of  the  rates  of  these  processes  shows  that  charge  exchange  is 
the  dominant  process,  at  least  for  the  species  that  have  an  ionization 
potential  low  enough  to  permit  exothermic  charge  exchange  with 
ionospheric  ions  (Caledonia,  et  al.,  1987). 

2.  Charge  exchange  processes  must  be  considered  separately,  even  if  they  are 
not  important  in  producing  the  seed  ionization.  There  are  two  different 
types:  charge  exchange  from  the  ambient  ions  (which  are  predominantly  0+) 
to  form  fast  ions 

O’*”  +  Mj  -♦  Mj'*’  +  0,  , 

which  has  a  cross  section  of  resonant  charge  exchange,  where 

one  beam  ion  charge  exchanges  with  a  fast  neutral  beam  molecule 

Ms+  +  Mf  Mf+  +  Mg,  ocer  » 

where  the  subscripts  refer  to  slow  and  fast  species.  This  process  is 
important  when  the  initial  ion  has  been  slowed  down  somewhat,  so  that  the 
resonant  charge  exchange  process  increases  the  directed  kinetic  energy 
along  the  beam  of  the  fast  ions  by  changing  a  slow  ion  into  a  fast  ion. 
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In  the  simulation,  the  ambient  charge  exchange  process  produces  one  new 
ion  from  the  beam  molecules  and  it  removes  one  ion  from  the  ambient  ion 
population.  The  charge  exchange  processes  do  not  change  the  number  of 
electrons. 

3.  Plasma  effects  will  heat  the  electrons  and  slow  the  fast  ions  by  develop¬ 
ing  an  electric  field  that  develops  as  a  result  of  exciting  a  plasma 
instability.  This  effect  is  discussed  more  in  other  sections,  and  this 
is  the  main  plasma  effect  that  we  wish  to  simulate.  The  heating  of  the 
electrons  is  essential  to  provide  the  mechanism  to  increase  the  charge 
density.  This  field  also  slows  the  beam  ions,  which  enhances  the 
resonant  charge  exchange  process. 

4.  Ionization  of  fast  neutrals  by  these  heated  electrons  is  the  mechanism  to 
increase  the  plasma  density.  This  process  provides  additional  fast  ions 
to  heat  more  electrons.  This  step  provides  feedback  that  can  produce  a 
discharge  under  conditions  where  each  fast  ion  transfers  more  energy  to 
the  plasma  electrons  than  the  average  energy  lost  by  the  electrons  to 
produce  another  ion  pair.  As  a  result  of  the  ionization,  there  are  now 
two  electrons  and  a  new  fast  ion. 

5.  Excitation  of  fast  neutrals  by  electron  collisions  decreases  the  energy 
of  the  plasma  electrons.  These  excited  species  can  be  ionized  by 
low-energy  electrons.  (Provided  that  the  electrons  have  an  energy 
greater  than  I  -  X,  where  I  is  the  ionization  potential  of  the  ground 
state  of  the  molecule  and  X  is  the  excitation  energy  of  the  excited  state 
formed.)  If  the  excited  species  dissociate  or  radiate  before  they  are 
ionized,  then  this  step  is  an  energy  sink  that  increases  the  average 
energy  lost  by  the  electrons  to  produce  an  ion  pair.  We  simulate 
conditions  where  the  excited  molecule  will  usually  radiate  before  it  is 
ionized,  so  our  simulation  considers  the  excitation  collisions  to  remove 
energy  from  the  electrons. 
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6. 


Elastic  collisions  between  the  electrons  and  the  fast  neutrals  will 
reorient  the  electron  velocity  vectors.  We  ignore  the  elastic  collisions 
of  the  electrons  with  the  ambient  species,  as  we  presume  that  the  beam 
number  density  is  much  larger  than  the  ambient  number  density  for  the 
regions  we  simulate.  When  the  neutral  density  is  large  enough,  there 
will  be  resistive  heating  of  the  electrons  (Machida  and  Goertz,  1986). 

In  our  simulations  the  neutral  number  densities  are  low  enough  that 
resistive  heating  is  not  important.  In  our  simulations  we  have  ignored 
the  small  amount  of  kinetic  energy  that  can  be  transferred  between  the 
massive  neutrals,  with  a  mass  m^j,  and  the  electrons,  with  mass  m^.  The 
fractional  energy  change  in  a  collision  is  of  the  order  of  (mg/mjij),  so 
the  major  effect  of  the  electron-neutral  elastic  collisions  is  to 
reorient  the  electron  velocity  vector. 

7.  Elastic  collisions  between  the  ions  and  the  fast  neutrals  will  reorient 
the  ion  velocity  vectors  and  transfer  kinetic  energy.  Again,  we  ignore 
the  elastic  collisions  with  the  ambient  neutral  species. 

3.5.2  Equations  for  CIV  Process 

In  this  section  we  give  equations  that  describe  the  CIV  process  in  the  pres¬ 
ence  of  the  molecular  physics  effects  of  charge  exchange,  ionization  by  elec¬ 
tron  collision,  electron  energy  loss  by  excitation  collisions,  and  elastic 
collisions  of  electrons  and  ions  with  neutrals.  We  consider  three  species  in 
the  plasma:  beam  ions  that  initially  have  a  directed  velocity  equal  to  that 
of  the  fast  neutral  beam,  ambient  ions  with  no  directed  velocity,  and  elec¬ 
trons.  Both  the  electrons  and  the  ambient  ions  have  thermal  energies.  We 
consider  that  the  neutral  gas  molecules  can  charge  exchange  with  the  ambient 
ions  to  form  beam  ions.  In  this  case,  we  can  neglect  other  ionization  proc¬ 
esses  to  form  the  initial  ions  (Caledonia,  et  al.,  1987).  We  start  with  an 
initial  fraction  of  the  total  number  of  ions  as  already  present  as  beam  ions. 

The  equation  for  the  density  of  the  beam  component  is 
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(101) 


dnb/dt  -  Vj  aj  ng  +  Vggg  rig  , 

where: 

vj  is  the  collision  frequency  for  ionizing  collisions 
aj  is  the  fraction  of  electrons  with  an  energy  greater  than  the  neutral 
molecule  ionization  potential  I,  and 
Vcea  Is  the  collision  frequency  for  charge  exchange  collisions  with  the 
ambient  ions. 

The  collision  frequencies  are  each  related  to  an  average  of  a  cross  section  ct 
times  the  electron  velocity  Vg 

V  =  <CT  Vg>  n[^  ,  (102) 

so  that  V  is  proportional  to  the  neutral  number  density. 

The  equation  for  the  ambient  ions  is 

dna/dt  =  -  Vggg  n^  .  (103) 

The  density  change  of  the  electrons  is  given  by 

dng/dt  =  vj  aj  ng  ,  (104) 

as  we  neglect  other  ionization  sources,  such  as  photoionization,  as  well  as 
loss  terms  such  as  recombination  and  diffusion.  We  maintain  charge  neutrality 

*^e  “  '^b  *^a  ' 

The  presence  of  energy  gains  from  charge  exchange,  energy  losses  from 
excitation  collisions  and  from  only  a  fraction  of  the  ion  energy  being 
transferred  to  the  electrons  will  change  the  critical  velocity.  Mobius  and 
coworkers  (1987)  derive  the  expression 
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Vb  >  Vc  [(1/ti  +  K)/(l  +  OjO-^  ,  (105) 

where  v^,  is  the  conventional  critical  velocity  (2I/mb)®'^,  tl  is  the  fraction 
of  the  fast  ion  kinetic  energy  that  is  transferred  to  the  electrons  (which  has 
a  value  of  0.8  for  their  derivation),  K  is  the  factor  for  excitation  losses, 
and  C  gives  the  charge  exchange  effects.  These  factors  are  given  by 

K  =  <Ojj  Vg>  /  <(Tj^  Vg>  ,  (106) 

where  and  crj  are  the  cross  sections  for  excitation  and  ionization,  respec¬ 
tively.  Note  that  we  should  include  as  excitations  only  those  excitations 
that  remove  energy  from  the  system  by  processes  such  as  radiation  and  dissoci¬ 
ation.  Excited  states  that  live  long  enough  to  be  ionized  are  not  energy 
losses,  as  these  states  have  lower  ionization  energies  (McNeil,  Lai,  and 
Murad,  1987).  The  charge  exchange  factor  is 

C  =  f-1-33  vb/<ffi  Ve>  ,  (107) 

where  f  =  nb/n^  is  the  fraction  of  positive  ions  that  are  in  the  beam. 

3.5.3  Description  of  Method  of  Computer  Simulation 

In  order  to  gain  insight  into  the  CIV  process,  we  performed  various  numerical 
simulations  of  plasmas.  The  simulations  use  several  approximations  to 
represent  the  physics  while  still  allowing  economical  calculations.  In  this 
section  we  discuss  the  approximations  that  apply  to  the  simulation  of  plasma 
effects  in  plasmas  where  atomic  and  molecular  physics  are  absent. 

3. 5. 3.1  Approximations  in  the  Plasma  Simulations.  These  simulations  used  a 
modification  of  the  ESI  code,  which  is  a  particle-in-cell  code  developed  by 
Birdsall  and  coworkers  at  the  University  of  California  (Birdsall  and  Langdon, 
1985;  Hasegawa  and  Birdsall,  1964).  The  original  ESI  code  considered  one 
dimension  in  position  space  and  two  dimensions  in  velocity  space  (e.g.,  x  for 
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position  and  Vj^  and  Vy).  As  discussed  earlier  we  have  modified  the  code  to 
consider  three  dimensions  in  velocity  space,  so  that  we  now  consider  four  of 
the  six  coordinates  in  phase  space.  The  modification  allows  us  to  vary  the 
angle  between  the  magnetic  field  and  the  velocity  vector  of  the  fast  ions. 
For  most  of  the  calculations  we  consider  the  velocity  vector  of  the  fast  ions 
to  be  parallel  with  the  electric  field  vector,  and  both  of  these  to  be 
parallel  to  the  one  positional  dimension  in  the  simulation. 

The  philosophy  used  by  Birdsall  is  that  the  simulations  are  intended  to 
produce  the  essence  of  the  plasmas,  without  all  of  the  details.  The  important 
physical  properties  that  the  simulations  have  to  preserve  and  the  approxima¬ 
tions  used  are  discussed  next. 

1.  The  plasmas  that  we  simulate  are  collisionless,  as  far  as  electron- 
electron  collisions  are  concerned.  This  can  be  be  seen  from  the  number 
of  particles  in  the  Debye  cube,  which  is  discussed  below. 

2.  We  are  interested  in  lengths  that  are  longer  than  the  Debye  length,  Xq 

Aq  =  [kT/(4rtnee2)|0-5 

=  743  T^'5  ng~®*5  cni 

where  T  is  the  electron  temperature  (eV)  and  n^  is  the  electron  number 
density  (1/cm^).  Thus,  for  T  =  0.1  eV  and  ng  =  10^/cm^,  Xq  =  0.23  cm. 

The  number  of  electrons  in  a  Debye  cube,  Np  is 

Nd  =  Ug  Xp3  . 

This  number  is  large,  1.2  x  lO'^,  for  the  example  with  ng  =  10^  cm^.  From 
Np  we  can  show  that  the  plasmas  are  collisionless,  as  far  as  electron- 
electron  collisions  are  concerned.  This  is  because  the  collisionless 
criterion  is 
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^e/“pe  :  ln(ND)/ND  «  1- 


where  Vg  is  the  electron-electron  collision  frequency  and  Wpg  is  the 
electron  plasma  frequency.  This  criterion  can  be  satisfactorily  met  for 
our  simulations  with  Nq  values  as  low  as  10.  Another  way  to  estimate  Vg 
is  to  use  the  Coulomb  cross  section  for  a  90  degree  deflection  (the  cross 
section  for  the  deflection  from  multiple  collisions  is  larger  than  the 
cross  section  for  the  deflection  in  a  single  collision),  which  is 

ogO  ~  2.6  X  10“^®/Eg  cm^, 

when  the  electron  kinetic  energy  Eg  is  given  in  eV  (Krall  and 
Trivelpiece,  1986,  p.  294).  Then 

\>e  =  ogQ  Vg  ng , 

which  is  1.5  to  150  collisions/s  for  Eg  =  1  eV  and  ng  =  10^  to  10^;  it  is 
a  factor  of  4  lower  when  Eg  =  16  eV.  Thus  the  assumption  of  collision¬ 
less  plasmas  seems  justified  for  the  simulations,  at  least  for  plasma 
densities  lower  than  10^,  although  it  may  not  be  justified  for  densities 
larger  than  10^.  This  is  because  the  growth  times  are  less  than  1  ms. 

3.  The  physical  picture  of  a  plasma  is  that  the  thermal  kinetic  energy  is 

large  in  comparison  with  the  potential  energy  associated  with  the  Coulomb 
forces  between  the  ions  and  the  electrons.  The  ratio 

Thermal  Kinetic  Energy/Microscopic  Potential  Energy  =  Nq. 

The  fundamental  physics  of  the  plasma  requires  that  this  ratio  is  much 
greater  than  one,  but  it  does  not  have  to  be  larger  than  10  for  the 
simulation  to  be  effective.  Thus  we  do  not  have  to  actually  use  a 
simulation  with  10^  electrons  in  a  Debye  cube;  we  can  get  by  with  only 
10.  The  plasma  parameter,  which  should  be  much  less  than  one,  is  the 
reciprocal  of  Np. 
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4.  In  the  simulations  we  use  a  spatial  grid  to  simplify  the  calculation  of 
the  fields.  This  has  the  effect  of  smoothing  away  information  on  the 
fields  and  forces  at  distances  less  than  one  grid  cell. 

5.  Periodic  boundary  conditions  effectively  make  the  length  of  the  system 
infinite.  This  correctly  embodies  the  physical  property  that  we  are 
interested  in  lengths  much  greater  than  the  Debye  length,  but  it  fails  to 
consider  effects  that  result  from  the  finite  size  of  the  system. 

6.  The  finite  size  of  the  ambient  plasma  and  of  the  fast  ion  beam  are 

ignored.  Also,  we  do  not  consider  the  density  variations  with  distance 

in  this  calculation. 

7.  We  assume  a  uniform  flow  of  neutral  particles  with  a  constant  flow 
velocity.  The  neutral  particle  number  densities  are  considered  to  be 
much  greater  than  the  ion  number  densities,  so  that  we  do  not  need  to 
follow  the  motion  of  the  individual  neutral  particles.  In  the 
simulations  we  always  consider  the  fast  neutrals  to  have  a  fixed  velocity 
of  8  km/s,  which  is  typical  of  the  velocity  of  the  Orbiter.  We  ignore 
the  perturbation  of  the  neutral  density  as  a  result  of  collisions  with 
neutral  particles  in  the  ambient. 

8.  The  ionization  of  a  fast  neutral  is  assumed  to  make  no  change  in  the 
momentum  of  the  heavy  particle,  so  that  the  fast  ion  has  the  same 
momentum  components  that  the  fast  neutral  had  before  the  ionization 
event.  This  is  a  good  approximation  for  ionization  by  electron  impact, 
and  it  is  still  a  reasonable  approximation  for  ionization  by  charge 
exchange. 

9.  The  simulations  use  artificially  small  ratios  of  the  mass  of  the  ion  to 
the  mass  of  the  electron,  mj/mg.  We  used  ratios  of  100  to  400,  rather 
than  ratios  of  5.1  x  10^  (for  N2‘'')  to  2.4  x  10^  (for  Xe"'').  This  speeds 
the  computation  while  still  preserving  the  essential  physics  that  the 


63 


ions  are  much  more  massive  than  the  electrons.  Also,  we  take  the  ions  to 
be  unmagnetized  (i.e.,  they  have  an  infinite  gyroradius),  which  correctly 
represents  another  important  physical  difference  between  the  ions  and  the 
electrons. 

10.  The  simulations  start  with  either  a  beam  of  fast  ions,  a  background  of 
thermal  ions,  and  a  background  of  thermal  electrons,  or  else  they  start 
without  the  background  of  thermal  ions.  In  all  cases  the  plasmas  are 
electrically  neutral.  That  is,  we  take  the  number  of  electrons  to  equal 
the  total  number  of  ions,  and  the  ions  are  singly  charged. 

11.  The  initial  velocity  distribution  of  the  background  electrons  and  of  the 
background  ions  are  each  assumed  to  be  Maxwellian  initially.  However, 
the  distributions  do  not  remain  Maxwellian  as  the  plasma  heats.  We 
follow  the  velocities  of  the  individual  particles. 

3. 5. 3. 2  Approximations  that  Should  be  Tested  in  Further  Work.  The  validity 
of  some  of  the  approximations  made  in  this  analysis  has  not  been  fully  tested. 
Their  limitations  should  be  evaluated  further  in  future  calculations.  This 
section  will  discuss  some  of  these  as  well  as  indicating  simple  modifications 
which  can  be  used  in  the  numerical  calculations  to  assess  their  significance. 

1.  Finite  geometry  and  spatial  variation  of  background  conditions  have  not 
been  considered.  Detailed  calculation  would  clearly  involve  a 
significantly  more  complex  computer  program.  For  a  finite  gas  cloud, 
however,  one  may  anticipate  that  both  wave  energy  and  heated  electrons 
may  leave  the  system.  If  this  happens  in  times  comparable  to  or  less 
than  the  times  required  for  the  waves  or  electrons  to  perform  their 
function,  converting  directed  energy  to  thermal  energy  and  then  ionizing 
neutra?s,  one  would  anticipate  that  the  process  would  become  less 
efficient. 
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The  effect  of  wave  energy  leaving  the  system  can  be  approximated  by 
assuming  a  damping  rate  equal  to  the  group  velocity  of  the  wave  divided 
by  the  size  of  the  region.  This  damping  rate  should  be  subtracted  from 
the  instability  growth  rate  and  can  be  incorporated  in  the  numerical 
analysis.  The  waves  of  interest  have  a  group  velocity  that  is 
predominantly  parallel  to  the  magnetic  field  and  is  proportional  to  the 
parallel  component  of  the  wave  vector.  For  typical  conditions  this 
damping  rate  is  not  negligible  for  waves  of  interest.  In  the  numerical 
calculations  the  angle  of  the  wave  vector  has  been  chosen  to  correspond 
to  the  maximum  growth  rate  in  an  infinite  uniform  medium.  Including 
damping  would  presumably  shift  the  maximum  growth  rate  to  smaller  values 
of  the  parallel  component  of  the  wave  vector  and  reduce  the  magnitude  of 
the  maximum  growth  rate.  The  overall  effect  on  rate  and  efficiency  of 
conversion  of  energy  to  ionization  should  be  calculated. 

The  effect  of  loss  of  high  energy  electrons  from  the  region  of  interest 
can  similarly  be  estimated  by  assuming  that  high  energy  electrons  leave 
the  region  at  a  rate  equal  to  the  electron  velocity  parallel  to  the 
magnetic  field  divided  by  the  size  of  the  region.  (At  high  neutral 
densities  the  electron  may  need  to  be  considered  to  diffuse  out  of  the 
region.)  This  loss  of  high  energy  electrons  which  are  then  replaced  by 
cold  ambient  electrons  can  be  added  to  the  computer  code  as  pseudo¬ 
collisions.  Rough  estimates  indicate  that  this  process  may  lead  to  a 
noticeable  reduction  in  the  efficiency  of  conversion  of  energy  to 
ionization. 

The  present  analysis  assumes  that  electron-electron  collisions  do  not 
occur.  This  is  certainly  valid  for  consideration  of  the  wave  growth  and 
conversi''n  of  energy  to  electron  kinetic  energy.  However,  since 
electrons  transfer  energy  very  efficiently  among  themselves,  the 
assumption  that  the  electron  distribution  function  becomes  non-Maxwellian 
is  not  valid  under  some  conditions  of  interest.  The  computer  code  can 
be  modified  to  assume  that  the  electrons  are  always  Maxwellian  but  that 
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the  electron  temperature  rises  as  the  electron  energy  increases.  The  net 
effect  of  this  on  ionization  efficiency  is  hard  to  predict  with 
precision.  The  non-Maxwellian  distribution  has  more  electrons  than  a 
Maxwellian  at  intermediate  energies  but  a  sharp  cutoff  at  high  energies. 
Since  the  ionization  rate  depends  on  the  number  of  electrons  above  the 
ionization  energy,  the  effect  on  the  rate  of  ionization  will  depend  on 
where  the  cutoff  is  compared  to  the  ionization  energy.  The  tail  of  the 
distribution  is  expected  to  extend  to  electron  energies  that  are  at  or 
above  the  ion  kinetic  energy  (Papadopolous,  1984),  which  is  above  the 
ionization  potential.  Therefore  the  Maxwellian  distribution  is  expected 
to  be  less  efficient  at  producing  ionization. 

3.5.4  Details  of  the  Simulation  of  Molecular  Physics 

In  this  section  we  shall  discuss  details  of  how  the  simulation  was  performed. 
The  object  of  the  simulations  was  to  learn  more  about  how  CIV  may  occur  in  a 
gas  release  from  Orbiter.  Thus,  it  is  important  to  use  realistic  cross 
sections  in  order  to  see  the  effects  that  occur  in  an  actual  gas.  For  this 
simulation,  we  used  data  from  nitrogen  molecules  to  represent  the  processes 
expected  for  different  gas  releases.  This  simulation  has  ignored  the  energy 
lost  to  vibrational  and  rotational  excitations.  These  losses  are  important 
when  the  electron  energy  is  less  than  the  energy  for  electronic  excitation. 
Thus  our  estimates  of  electron  heating  will  exceed  that  in  a  real  molecular 
gas.  The  overestimate  is  less  important  for  the  low-pressure,  high  neutral 
beam  ene-ry  runs  where  a  significant  fraction  of  the  electrons  are  quickly 
heated  above  the  ionization  energy.  We  now  give  details  of  the  cross  section 
and  angular  distribution  calculations. 

3. 5. 4.1  Method  of  Simulating  Electron-Neutral  Collisions.  The  procedure  is 
to  run  the  plasma  simulation  code  and  allow  the  electrons  to  heat.  At 
periodic  intervals  we  call  a  subroutine  to  determine  if  an  electron  collision 
(or  other  molecular  physics)  has  occurred.  This  is  done  by  using  the  total 
cross  section  for  electron-neutral  collisions, 


66 


•^tot  =  ‘^inel  +  ®el  * 


which  is  the  sum  of  the  inelastic  and  elastic  cross  sections  calculated  for  an 
electron  with  an  energy  Eg.  This  total  cross  section  is  used  with  the 
electron  velocity  Vg  (but  scaled  to  the  velocity  of  an  actual  electron,  rather 
than  the  velocity  of  the  simulation  electron),  and  the  number  density  of  the 
fast  neutral  beam  nj^j  to  compute  the  mean  free  time  between  collisions 

Tc  =  l/(fftot*^N'^e^ 

This  Tg  value  determines  the  distribution  of  times  between  collisions.  Ve  can 
sample  this  distribution  by  finding  a  random  number  R  selected  from  a  uniform 
distribution  from  0  to  1  and  using  the  equation 

At{^  =  -  ln(R)/Tg 

If  AtR  is  greater  than  the  time  that  has  elapsed  since  the  last  subroutine 
call,  then  the  simulation  considers  that  the  electron  in  question  had  a 
collision.  At  each  collision,  another  random  number  is  selected  to  determine 
the  nature  of  the  collision.  If  this  number  is  less  than  the  fraction  of 
collisions  that  give  elastic  collisions  (which  is  equal  to  ,  then  the 

collision  is  elastic.  If  the  number  is  larger,  then  the  collision  is 
inelastic  and  the  random  number  is  examined  to  see  which  inelastic  process  is 
excited.  Here  is  considered  to  be  given  by 

^tot  ~  ^el  ^  » 

where  oj  (Eg,Ej)  is  the  inelastic  cross  section  for  transfer  of  an  energy  Ej 
to  the  neutral  molecule.  In  practice  we  only  used  a  limited  number  of  Ej 
values,  with  one  value  (12  eV)  to  simulate  excitation  and  several  values  to 
represent  ionization.  These  values  were  taken  to  be  the  ionization  potential 
(I  =  15.6  eV)  plus  secondary  electron  energies  of  1,  2,  4,  7,  11,  16,  22,  30, 
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AO,  60,  and  80  eV.  The  maximum  secondary  electron  energy  in  a  collision  of  an 
electron  with  an  initial  energy  of  Eg  is  given  by  0.5(Eg  -  I). 

If  the  collision  is  elastic,  we  find  the  new  direction  cosines  of  the  electron 
by  taking  two  more  random  numbers,  one  to  give  the  azimuthal  scattering  angle 
(which  is  uniformally  distributed  between  0  and  2ii),  and  one  to  give  the 
deviation  from  the  original  direction  0.  The  0  angle  is  found  by  integrating 
the  angular  distribution  (Ee>9)  times  sin(0)d0  to  give  a  uniform 
probability  between  0  and  the  total  elastic  cross  section.  The  integrated 
distribution  is  used  with  a  random  number  to  select  the  value  of  the 
scattering  angle  0.  In  this  way  we  allow  backseat tering  to  occur,  as  0  varies 
between  0  and  n. 

If  the  collision  is  inelastic,  we  find  the  new  direction  cosines  in  a  similar 
manner,  except  we  use  the  inelastic  angular  distribution  aj  (Eg,Ej)  for 
scattering  by  an  inelastic  collision  that  transfers  an  energy  E j . 

The  cross  sections  are  functions  of  the  kinetic  energy  of  the  electron. 
Actually,  they  depend  on  the  kinetic  energy  of  the  relative  motion  of  the 
electron  and  the  neutral  molecule.  The  electrons  have  much  greater  velocities 
than  the  molecules,  even  for  molecules  traveling  at  8  km/s,  as  the  electrons 
travel  at  593  Eg^*^  km/s,  where  Eg  is  the  electron  kinetic  energy  in  eV.  Thus 
an  electron  at  20  eV  is  traveling  at  2.7  x  10^  km/s.  In  addition,  the  actual 
mass  of  the  ion  is  5  x  10^  (or  more)  times  as  great  as  the  mass  of  the 
electron.  Therefore,  we  need  only  consider  the  motion  of  the  electron.  In 
the  simulation,  we  use  cross  sections  based  on  actual  data  to  give  a  realistic 
energy  dependence. 

The  physical  inputs  into  the  simulation  are  the  values  of  absolute  cross 
sections  ffgi  (E)  and  ffinelj  (Ee,Ej),  the  relative  cross  sections  for  angular 
scattering  Ogj  (Eg,0)  and  ffinelj  (Ee*Ej,6),  the  relative  cross  section  for 
producing  a  fast  secondary  electron  at  an  angle  0^,  Osec  (Ee>Ej,6s)’ 
next  subsections,  we  shall  discuss  our  choices  for  each  of  these  values. 
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3. 5. 4. 2  Absolute  Cross  Sections  for  Elastic  and  Inelastic  Scattering.  The 

absolute  cross  sections  (together  with  the  gas  number  density)  are  the 
parameters  that  determine  the  length  of  the  path  between  collisions  as  well  as 
the  distribution  of  energy  losses  produced  in  inelastic  collisions.  We  need  a 
general  prescription  for  determining  agi  (E^)  and  ffjnelj  wide 

ranges  of  Eg  and  E j .  For  this  purpose  we  have  used  the  expression  used  by 
Jackman  and  Green  (1979)  (JG)  for  Ogj^  and  the  expression  given  by  Porter, 
Jackman,  and  Green  (1976)  (PJG)  for  Oinel*  These  expressions  give  values  that 
compare  with  the  recommended  values  given  by  Itikawa,  et  al.  (1986)  to  within 
10  to  20  percent.  We  have  also  compared  the  differential  cross  section, 
doionj  producing  a  secondary  electron  with  an  energy 

(Ej-15.6)  eV  with  Fig.  9.8  in  Itikawa,  et  al.  (1986).  Here  our  values  usually 
are  some  5  to  15  percent  high. 

3. 5. 4. 3  Elastic  Scattering  Angular  Distribution.  The  relative  differential 

scattering  cross  section  Ogi  (Eg,  0)  determines  the  distribution  of  the 
angular  scattering  by  elastic  processes  where  the  election  energy  Eg  is 
unchanged  during  the  collision.  We  have  used  either  the  general  prescription 
for  Ogj^  (Eg,0)  given  by  Jackman  and  Green  (1979)  or  a  modification  of  this 
distribution.  Some  simulations  were  run  with  a  modification  where  the 
contribution  from  small  angles  is  diminished.  The  modification  was  to  make 
®el  constant  for  0  <  3  deg.  The  simulations  reported  here  take  the 

angular  distribution  to  be  the  same  for  all  electron  energies  below  17  eV,  and 
they  use  the  distribution  calculated  for  17  eV.  The  energy  dependence  of  the 
angular  distribution  is  considered  for  electron  energies  greater  than  17  eV. 
The  neglect  of  the  energy  dependence  below  17  eV  will  reduce  the  effects  of 
elastic  collisions,  as  the  true  angular  distribution  will  have  a  smaller 
contribution  from  the  forward  peak. 

3. 5. 4. 4  Inelastic  Scattering  Angular  Distribution.  The  inelastic  scattering 
can  also  be  simulated  by  knowing  the  absolute  total  cross  section  and  the 
relative  angular  distribution.  Here  we  focus  on  the  relative  values  for 
Oineij  (Ez»Ej,0).  The  inelastic  scattering  of  electrons  can  be  described  by 
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using  the  generalized  oscillator  strength  (GOS)  (Inokuti,  1971).  A  complete 
knowledge  of  of  the  GOS  provides  a  method  of  defining  the  absolute  angular 
scattering  cross  section.  The  GOS  has  useful  properties,  such  as  going  to  the 
optical  oscillator  strength  in  the  limit  of  zero  momentum  transfer,  and  there 
are  general  sum  rules  for  various  moments  (Inokuti,  1971).  Porter,  Jackman, 
and  Green  (PJG)  considered  these  properties  in  order  to  fit  a  parameterized 
form  of  the  Massey-Mohr-Bethe  surface  (which  is  based  on  the  theoretical  form 
of  GOS  that  is  known  to  be  correct  for  H  atoms).  The  fit  is  given  in  their 
equation  (14)  and  in  their  Table  II,  with  the  following  corrections  (Green, 
personal  communication,  1985);  a  constant  of  0.654  is  added  to  82  ® 

constant  of  1.0  is  added  to  y  .  Our  simulations  use  the  total  inelastic 
a(E,Eg)  given  by  PJG  together  with  the  relative  angular  distributions  using 
a(E,Eg)  derived  from  their  GOS  surface. 

3. 5. 4. 5  Method  of  Simulating  Ion-Neutral  Collisions.  The  ion-neutral 
collisions  are  simulated  in  a  similar  manner  to  the  electron-neutral 
collisions.  In  both  cases,  we  assume  that  the  neutral  number  density  is 
uniform,  constant,  and  much  larger  than  the  density  of  ambient  species.  For 
this  case  the  subroutine  call  determines  if  an  ion  collision  has  occurred. 

This  is  done  by  using  the  total  cross  section  for  ion-neutral  collisions, 

'^ti  =  ^^ce  '^eli» 

which  is  the  sum  of  the  charge  exchange  and  elastic  cross  sections  calculated 
for  an  ion  with  an  energy  Ej .  The  cross  section  for  charge  exchange  with 
ambient  O'*"  ions  is  estimated  from  the  rate  for  this  reaction  for  a  relative 
velocity  of  8  km/s  (Albritton,  1978)  as  5  x  10"^^  cm^.  The  cross  section  for 
resonant  charge  transfer  is  needed  as  a  function  of  relative  velocity;  we 
approximated  the  cross  sections  from  Stebbings  and  coworkers  (1963)  as  a 
linear  function  of  the  velocity,  with  Ocer  =  33  x  10~^^  cm^  at  v,.  =  8  km/s  and 
the  slope  =  -2.4  x  10“^^  cm2/(cm/s).  The  elastic  cross  sections  were 
calculated  from  an  evaluation  (Mahan  and  Person,  1964)  corresponding  to 
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scattering  from  collisions  where  the  potential  energy  P  is  represented  as  the 
sum  of  an  ion-induced  dipole  attraction  and  a  hard-sphere  core 

P  =  ap  e2  /(2  R^)  R  >  s 

P  =  “  R  <  s 

The  hard  sphere  radius  s  was  estimated  to  be  2.55  x  10~®  and  3.0  x  10”®  cm  for 
O'*"  and  N2''',  respectively,  by  analogy  with  Goldschmidt's  values  (Loeb,  1960). 
The  polarizability  ap  of  the  neutral  is  approximated  as  1.76  x  10”^'^  cm^,  the 
value  for  N2. 

The  total  cross  section  is  used  with  the  ion  velocity  vj  (which  was  scaled  to 
the  velocity  of  an  actual  O’*"  or  N2‘*^  ion),  and  the  number  density  of  the  fast 
neutral  beam  nj,}  to  compute  the  distribution  of  times  between  collisions.  As 
for  the  electrons,  we  use  a  random  number  to  determine  whether  the  ion  in 
question  had  a  collision.  At  each  collision,  another  random  number  is 
selected  to  determine  the  nature  of  the  collision.  If  this  number  is  less 
than  the  fraction  of  collisions  that  give  elastic  collisions  (which  is  equal 
to  ffeli^'^ti)’  collision  is  elastic.  If  the  number  is  larger,  then 

the  collision  gives  charge  exchange,  which  is  either  from  the  ambient  ion  or 
resonant,  depending  on  which  type  of  ion  we  are  following. 

If  the  collision  is  elastic,  the  effect  of  the  elastic  collision  is  to  rotate 
the  relative  velocity  vector  through  some  angle  X  in  the  center  of  mass 
system.  The  distribution  of  X  favors  small  values  of  X>  and  random  choices  of 
X  used  the  angular  distribution  calculated  from  the  tabulated  results  for  the 
Langevin  potential  given  above  (Mahan  and  Person,  1964).  The  choice  of  X 
determines  the  kinetic  energy  transferred  in  the  collision.  In  addition  we 
take  another  random  number  to  give  the  azimuthal  scattering  angle  which 
then  determines  the  direction  cosines  of  the  new  ion.  Additional  details  on 
the  elastic  scattering  calculation  can  be  found  in  the  report  by  Caledonia, 
Person,  and  Hastings  (1986).  If  the  collision  is  charge-exchange,  we  find  the 
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new  direction  cosines  by  assuming  that  the  new  fast  ion  has  the  same  velocity 
as  that  of  the  fast  neutral  beam. 


3.5.5  Results  of  Computer  Simulations 

The  simulation  uses  enough  molecular  physics  to  allow  study  of  several  of  the 
competing  effects  that  are  related  to  the  CIV  process.  The  number  density 
(pressure)  of  the  neutral  beam  and  the  kinetic  energy  of  the  fast  ions  were 
the  major  parameters  that  we  varied.  There  are  competing  effects  that  result 
from  increasing  the  pressure  of  the  neutrals.  One  effect  is  to  increase  the 
number  of  elastic  collisions  between  electrons  and  neutrals.  These  collisions 
scatter  the  electrons,  which  reorient  their  velocity  vectors.  The  effect  of 
these  collisions  is  to  reduce  the  electron  heating  via  the  interaction  with 
the  electric  field.  In  these  simulations  the  effect  of  the  elastic  scattering 
of  the  electrons  is  reduced  because  the  low  energy  electrons  have  an  angular 
distribution  with  too  large  a  peak  in  the  forward  direction.  Approximately 
AO  percent  of  the  elastic  collisions  result  in  deflections  that  are  less  than 
30  deg.  In  addition,  the  neglect  of  vibrational  and  rotational  excitation 
collisions  will  also  increase  the  electron  heating  in  these  simulations. 

Another  effect  of  increasing  the  neutral  number  density  is  to  increase  the 
rate  of  charge  exchange,  both  with  the  ambient  ions  and  with  the  ions  in  the 
beam.  The  resonant  charge  exchange  with  ions  in  the  beam  has  a  larger  cross 
section  in  the  simulation  (~  33  x  10~^^  versus  5  x  10“^^  for  ambient  ions), 
but  the  relative  velocity  for  this  process  is  essentially  zero  unless  the  beam 
ion  has  been  slowed  somewhat.  Thus  the  resonant  charge  exchange  process  is 
more  likely  to  occur  after  the  electric  field  has  built  up  and  some  of  the 
fast  beam  ions  have  lost  kinetic  energy.  The  resonant  charge-exchange 
process  will  increase  the  energy  of  the  beam  ions.  Charge  exchange  with  the 
ambient  ions  will  increase  the  fraction  of  the  ions  that  are  in  the  fast  ion 
beam.  Increasing  this  fraction  increases  the  electron  heating  via  plasma 
effects,  which  increases  the  ionization  via  hot  electron  collisions.  The 
change  of  this  fraction  will  change  the  k^jgj^  value  for  the  fastest  growing 
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waves  of  the  instability.  Changing  kn,^^  reduces  the  rate  of  growth  of  the 
electric  field  energy  in  comparison  to  when  the  fraction  of  fast  ions  is 
constant.  On  the  other  hand,  when  the  fraction  of  beam  ions  is  constant  there 
are  other  resonance  effects  that  cause  the  electric  field  energy  to  go  through 
a  maximum.  This  is  shown  by  the  upper  curve  in  Fig.  16,  which  is  for  a  run 

where  the  neutral  density  is  1000  times  smaller  than  the  lower  curve.  The 

lower  curve  shows  the  electric  field  energy  for  a  case  where  charge  exchange 
occurs,  which  keeps  'hanging  the  fraction  of  beam  ions.  Here  we  see  that  the 
maximum  field  energy  is  less.  However,  the  field  energy  does  not  have  a 
maximum.  The  typical  behavior  is  for  the  field  energy  to  reach  a  plateau 
value  (which  sometimes  occurs  rather  early),  or  else  the  field  continues  to 
slowly  increase  after  an  initially  fast  build-up. 

Table  1  gives  a  summary  of  results  from  the  simulations.  The  notes  to  the 
table  contain  details  about  the  simulations.  In  particular  note  that  the 
times  are  given  in  an  internal  time  unit.  For  a  magnetic  field  of  0.45  Gauss, 
one  unit  is  2.16  us  for  O'*’,  2.85  us  for  N2^,  2.95  us  for  NO"*',  3.58  us  for 

C02'*',  and  6.17  us  for  Xe"*".  As  a  result  of  this  time  unit,  the  number 

densities  (l/cm^)  of  beam  neutral  molecules  are  multiplied  by  a  mass-dependent 
factor  Fp.  The  value  of  Fp  is  84.3  for  O'*",  112  for  N2"'',  116  for  NO"*^,  140  for 
C02^,  and  241  for  Xe"*".  In  this  discussion  we  shall  consider  thr>*  Fp  is 
approximately  100,  as  the  value  of  Fp  is  84  to  113  for  ions  around  the  masses 
used  in  the  simulation,  and  we  take  the  time  unit  as  approximately  2.5  us  so 
that  the  times  of  461  and  692  correspond  to  1.2  and  1.7  ms,  respectively. 

Note  also  that  the  runs  initially  have  1024  electrons  and  1024  positive  ions. 
The  table  is  sorted  by  Egjj/I,  nj^,  and  f|j  (t=0).  The  simulations  all  assume 
Vjj  =  8  km/s,  so  the  change  in  is  the  result  of  changing  the  ionization 

potential  and  mass  of  the  beam  ion.  The  column  headings  are  explained  at  the 
end  of  the  table. 

3. 5. 5.1  Fraction  of  Electrons  Heated.  Table  1  shows  that  the  maximum 
fraction  of  electrons  heated  to  an  energy  greater  than  1/2,  fg]^,  and  the 
maximum  fraction  with  an  energy  greater  than  I,  fg2,  are  both  functions  of  the 
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ENERGY  (INTERNAL  UNITS) 


TABLE  1.  Summary  of  simulations  for  fast  ion  beam  traversing  background 
plasma. 
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TABLE  1.  Summary  of  simulations  for  fast  ion  beam  traversing  background 
plasma  (Continued). 
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TABLE  1.  Summary  of  simulations  for  fast  ion  beam  traversing  background 
plasma  (Continued). 
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TABLE  1.  Summary  of  simulations  for  fast  ion  beam  traversing  background 
plasma  (Continued). 


Notes : 

1.  All  runs  used  a  mass  ratio  of  the  positive  ion  to  electron  of  100. 

2.  The  plasma  simulation  considers  the  beam  ions  and  the  ambient  ions  to  have 
the  same  mass,  but  the  ion-neutral  collisions  use  a  mass  corresponding  to 
O'*’  for  the  ambient  ions  and  N2^  for  the  beam  ions. 

3.  The  beam  molecules  all  have  a  velocity  of  8  km/s. 

4.  The  initial  thermal  energy  of  electrons,  ambient  ions,  and  beam  ions  is 
0.1,  0.1,  and  0.01  eV,  respectively. 

5.  If  there  is  a  zero  entry  in  the  column  for  ion-elastic  collisions,  or  for 
charge  exchange,  that  is  an  indication  that  the  corresponding  cross 
section  has  been  set  to  zero. 

6.  All  runs  used  256  grid  cells  with  a  length  chosen  to  maximize  growth  for 
mode  one  in  the  linear  dispersion  analysis.  For  the  runs  where  fbO  is 
0.005  and  where  there  is  charge  exchange  with  ambient  ions,  the  linear 
anlysis  used  f  =  0.05  to  allow  for  the  growth  of  fb.  Run  60  is  a  repeat 
of  run  58  using  the  length  calculated  from  fbO  =  0.5;  this  shows  effects 
of  different  lengths. 

7.  All  runs  have  the  electron  plasma  frequency  2.26  times  the  electron 
gyrof requency ,  which  corresponds  to  an  electron  number  density  of 
10^  electrons/cm^  and  a  magnetic  field  of  0.A5  Gauss. 

8.  All  times  are  given  in  an  internal  time  unit.  For  a  magnetic  field  of 
0.45  Gauss,  one  unit  is  2.16  ys  for  O"*",  2.85  ys  for  N2‘'^,  2.95  ys  for  N0+, 
3.58  ys  for  C02''',  and  6.17  ys  for  Xe+.  As  a  result  of  this  time  unit,  the 
number  densities  (cm“3)  of  beam  neutral  molecules  are  multiplied  by  a 
mass-dependent  factor  Fp.  The  value  of  Fp  is  84.3  for  O'*',  112  for  N2''', 

116  for  N0+,  140  for  C02-^,  and  241  for  Xe+, 

9.  All  runs  initially  have  1024  electrons  and  1024  positive  ions,  with 
the  positive  ions  distributed  between  beam  ions  and  ambient  ions. 

10.  Runs  numbered  35  and  smaller  used  a  flattened  angular  distribution  for 
electron-neutral  elastic  collisions.  Runs  36  and  37  are  repeats  of 
runs  35  and  34. 
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TABLE  1.  Summary  of  simulations  for  fast  ion  beam  traversing  background 
plasma  (Concluded). 


11.  Runs  numbered  20  and  smaller  used  seven  energy  bins  for  the  electron 
energies,  rather  than  11  bins  as  used  for  later  runs.  Run  22  is  a 
repeat  of  run  15  to  show  the  minor  effects  of  this  change. 

Explanation  of  columns: 

Run  -  Run  number. 

Ekb/I  -  Kinetic  energy  of  beam  ion  divided  by  ionization  potential  I  of  beam 
molecule. 

nN*Fp  -  Number  density  (cm~^)  of  beam  neutral  molecules  multiplied  by  a 
mass-dependent  factor  Fp.  At  0.45  Gauss,  the  value  of  Fp  is 
84.3  for  0+,  112  for  N2+,  116  for  N0+,  140  for  CO2+,  and  241  for  Xe+. 
fbO  -  Fraction  of  ions  that  are  initially  in  the  beam,  fbO  =  n^/Og. 

tli  -  Time  before  first  one  of  electrons  has  more  than  1/2  of  kinetic 

energy.  See  note  8. 

tlf  -  Time  before  at  least  one  of  electrons  always  has  more  than  1/2  of 
kinetic  energy.  See  note  8. 

t2i  -  Time  before  first  one  of  electrons  has  more  than  I  of  kinetic  energy. 
See  note  8. 

t2f  -  Time  before  at  least  one  of  electrons  always  has  more  than  I  of 
kinetic  energy.  See  note  8. 

tion  -  Time  before  first  electron  ionization.  See  note  8. 
felmx  -  Maximum  fraction  (in  this  time  interval)  of  electrons  with  more  than 
1/2  of  kinetic  energy. 

fe2mx  -  Maximum  fraction  (in  this  time  interval)  of  electrons  with  more  than 
I  of  kinetic  energy, 
t  -  Time,  for  results.  See  note  8. 

Nex  -  Number  of  electronic  excitation  events.  See  note  9. 

Nion  -  Number  of  ionization  events  by  electron  collision.  See  note  9. 

Nele  -  (1/1000)  times  number  of  elastic  electron-neutral  collisions.  See 
note  9. 

Nelia  -  (1/1000)  times  number  of  elastic  ambient  ion-neutral  collisions.  See 
note  9. 

Nelib  -  (1/1000)  times  number  of  elastic  beam  ion-neutral  collisions.  See 
note  9. 

Ncea  -  Number  of  charge  exchange  events  with  ambient  ions.  See  note  9. 

Ncer  -  Number  of  resonant  charge  exchange  events  (with  slow  beam  ions).  See 

note  9. 

fEf  -  Fraction  of  total  energy  in  the  electric  field. 

fEke  -  Fraction  of  total  energy  in  the  kinetic  energy  of  the  electrons. 

fEkb  -  Fraction  of  total  energy  in  the  kinetic  energy  of  the  beam  ions. 

fEka  -  Fraction  of  total  energy  in  the  kinetic  energy  of  the  ambient  ions. 
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neutral  number  density  nj^,  the  initial  kinetic  energy  of  a  beam  ion  and 

the  initial  fraction  of  ions  in  the  beam  f^o-  Figures  17  and  18  are  plots  of 
fgl  and  fg2  versus  njij  for  times  around  1.2  ms.  The  plots  are  line  segments 
connecting  the  points  for  the  same  and  fjjQ  values.  The  lines  for  fjjQ  = 

0.5  show  that  more  electrons  are  heated  when  Egj,  is  twice  I  than  when  Egjj/I  is 
0.6  or  1.  Also,  increasing  nfij  reduces  the  fraction  heated,  especially  for 
Fig.  18,  and  there  seems  to  be  a  minimum  around  1  to  3  x  10^^  cm~^.  For  the 
Egb/I  =  2  case,  nj^  has  less  effect,  but  there  still  is  a  falloff  of  the  hotter 
electrons  in  Fig.  18  as  n^j  increases.  When  we  start  with  a  weakly  ionized 
beam  with  fj^g  =  0.005,  we  find  no  heating  until  there  are  enough  charge 
exchange  collisions  to  increase  the  number  of  beam  ions  and  to  pump  energy 
into  the  beam  via  resonant  charge  exchange  collisions.  Thus  Figs.  17  and  18 
show  these  cases  to  rise  steeply  with  increasing  njyj.  They  also  show  that  the 
fraction  of  electrons  heated  increases  as  Egj,  increases. 

3. 5. 5. 2  Rate  of  Electron  Heating.  The  results  on  the  fraction  of  electrons 

heated  can  be  understood  as  the  result  of  the  different  rates  of  heating  the 

electrons  for  different  neutral  number  density  nj^,  kinetic  energy  of  the  beam 

Egjj,  and  initial  fraction  of  ions  in  the  beam  f^g.  Figures  19  and  20  are 
plots  of  the  reciprocal  of  the  time  (in  internal  units)  to  heat  the  first 
0.001  of  the  electrons  to  a  kinetic  energy  of  1/2  or  I,  respectively.  The 
plots  are  line  segments  connecting  the  points  for  the  same  Egjj  and  f5g  values. 

The  lines  for  f^g  =0.5  show  that  more  electrons  are  heated  faster  when  Eg^^I 

is  two  than  when  Egjj/I  is  0.6.  Also,  increasing  nj^j  has  only  a  modest  effect 
when  f^g  is  0.5;  it  decreases  the  heating  rate  to  reach  an  energy  of  I  for  the 
Egb/I  =  0.6  case.  When  we  start  with  a  weakly  ionized  beam  with  fjjg  =  0.005, 
we  find  that  the  heating  rate  increases  with  n(^J,  which  increases  the  rate  of 
charge  exchange  reactions.  Figures  19  and  20  show  also  that  the  heating  rate 
increases  as  Egjj  increases.  Note  that  the  values  at  0.00204  and  0.00145 
represent  upper  li  its. 

3. 5. 5. 3  Disposition  of  the  Energy  of  the  Fast  Ion  Beam.  If  we  consider  the 
case  where  we  have  an  initial  fraction  of  ions  in  the  fast  ion  beam  and  we  do 
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Figure  17.  Maximum  fraction  of  electrons  heated  to  an  energy  of  1/2  during  the 
the  first  461  time  units  as  a  function  of  the  number  density  of 
neutral  molecules.  The  lines  are  labeled  with  the  energy  of  the 
beam  ions  Egj,  and  with  the  initial  fraction  of  ions  in  the  beam 
f{,Q.  The  solid  line  labeled  NO-IEL  represents  runs  where  there 
are  no  elastic  collisions  of  the  ions  and  neutrals.  Note  8  to 
Table  1  explains  the  mass-dependent  values  of  the  time  unit, 
which  is  2.85  us  for  N2^,  and  Fp,  which  is  112  for  N2'''  when  nj^j 
is  in  cm"^. 


Figure  18.  Maximum  fraction  of  electrons  heated  to  an  energy  of  I  during  the 
the  first  461  time  units  as  a  function  of  the  number  density  of 
neutral  molecules.  The  lines  are  labeled  with  the  energy  of  the 
beam  ions  Egj,  and  with  the  initial  fraction  of  ions  in  the  beam 
ffjQ.  The  solid  line  labeled  NO-IEL  represents  runs  where  there 
are  no  elastic  collisions  of  the  ions  and  neutrals.  Note  8  to 
Table  1  explains  the  mass-dependent  values  of  the  time  unit, 
which  is  2.85  ps  for  N2^,  and  Fp,  which  is  112  for  N2‘*^  when  njij 
is  in  cm"3. 
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1/TIME  (E  >  1/2) 


Figure  19.  Rate  of  heating  0.001  of  the  electrons  to  an  energy  of  1/2  as  a 
function  of  the  number  density  of  neutral  molecules.  The  lines 
are  labeled  with  the  energy  of  the  beam  ions  Egjj  and  with  the 
initial  fraction  of  ions  and  neutrals.  Note  8  to  Table  1 
explains  in  the  beam  f{jQ.  The  solid  line  labeled  NO-IEL 
represent.s  runs  where  there  are  no  elastic  collisions  of  the 
ions  and  neutrals.  Note  8  to  Table  1  explains  the  time  unit, 
which  is  2.85  ys  for  N2^,  and  Fp,  which  is  112  for  N2'''  when  n[,j 
is  in  cm“3. 
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1/TIME  (E  >  I) 


Figure  20.  Rate  of  heating  0.001  of  the  electrons  to  an  energy  of  I  as  a 

function  of  the  number  density  of  neutral  molecules.  The  lines 
are  labeled  with  the  energy  of  the  beam  ions  and  with  the 
initial  fraction  of  ions  in  the  beam  fjjQ.  The  solid  line 
labeled  NO-IEL  represnts  runs  where  there  are  no  elastic 
collisions  of  the  ions  and  neutrals.  Note  8  to  Table  1  explains 
the  mass-dependent  values  of  the  time  unit,  which  is  2.85  ps 
for  N2‘*^,  and  Fp,  which  is  112  for  N2'*^  when  njj  is  in  cm~^. 
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not  allow  charge  exchange  or  ion-neutral  elastic  collisions,  then  we  can 
determine  the  disposition  of  the  energy  that  was  initially  present  as  kinetic 
energy  of  the  fast  ion  beam.  This  energy  is  divided  between  the  field  energy 
and  the  kinetic  energy  of  (1)  the  electrons  Eg,  (2)  the  ambient  ions  Egg,  and 
(3)  the  beam  ions  Egij.  We  choose  the  angle  between  the  magnetic  field  and  the 
velocity  vector  of  the  beam  ions  to  maximize  the  electron  heating  in  the 
linear  dispersion  analysis.  At  this  angle  approximately  52  to  60  percent  of 
the  initial  kinetic  energy  of  the  beam  ions  is  transferred  to  Eg  and  4  to 
12  percent  to  Eg.  Thus  h  is  0.52  to  0.6  in  this  simulation,  lower  than  the 
0.8  estimated  by  Mobius  and  coworkers  (1987),  but  consistent  with  the  limit  of 
2/3  given  by  Galeev  (1981)  and  with  the  results  of  other  simulations  (e.g., 
McBride,  et  al,  1972;  Machida  and  Goertz,  1986). 

Figures  21  through  23  shows  the  fraction  of  the  total  energy  in  the  kinetic 
energy  of  the  electrons,  the  beam  ions,  and  the  ambient  ions,  respectively, 
all  evaluated  at  a  time  near  1.2  ms  (461  internal  units).  At  low  pressures 
the  fraction  in  electron  kinetic  energy  is  approximately  independent  of  nj^j. 
Here  it  is  given  by  0-  At  higher  pressures  energy  is  transferred  to  the  beam 
ions  faster  than  the  ions  transfer  energy  to  the  electrons.  Here,  the 
fraction  of  energy  in  the  beam  ions  increases  and  the  fraction  in  the 
electrons  decreases,  as  shown  in  Figs.  21  and  22.  Note  that  the  larger  values 
of  Egjj/I  give  a  larger  fraction  in  the  energy  of  the  beam.  This  is  the  result 
of  the  resonant  charge  exchange  energy  feeding  energy  into  the  beam  faster 
than  the  energy  is  transferred  to  the  electrons.  That  is,  the  larger  Egjj/I 
beams  heat  electrons  faster  (as  shown  in  Figs.  19  and  20).  As  a  result  the 
beam  ions  are  slowed  faster,  which  gives  more  resonant  charge  exchange 
collisions.  The  figures  are  not  adjusted  for  the  changing  proportion  of  the 
ions  that  are  in  the  beam.  This  explains  why  the  f^o  =0.5  cases  are  high  for 
their  value  of  Ef(^jj/I.  Also,  this  explains  the  decrease  in  the  ambient  ion 
energy  fraction  as  the  number  of  ambient  ions  decreases  for  n^yj  greater  than 
3  X  10^^  cm“3. 
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Figure  21.  Fraction  of  the  total  energy  in  the  electron  kinetic  energy 
after  461  time  units  as  a  function  of  the  number  density  of 
neutral  molecules.  The  lines  are  labeled  with  the  initial 
energy  of  the  beam  ions  and  with  the  initial  fraction  of 
ions  in  the  beam  ff,Q.  The  solid  line  labeled  NO-IEL 
represents  runs  where  there  are  no  elastic  collisions  of  the 
ions  and  neutrals.  Note  8  to  Table  1  explains  the 
mass-dependent  values  of  the  time  unit,  which  is  2.85  ys  for 
N2^,  and  Fp,  which  is  112  for  N2^  when  njq  is  in  cm“^. 


BEAM  ION  KE/TOTAL 


Figure  22.  Fraction  of  the  total  energy  in  the  beam  ion  kinetic  energy 
after  461  time  units  as  a  function  of  the  number  density  of 
neutral  molecules.  The  lines  are  labeled  with  the  initial 
energy  of  the  beam  ions  and  with  the  initial  fraction  of 
ions  in  the  beam  f^Q.  The  solid  line  labeled  NO-IEL 
represents  runs  where  there  are  no  elastic  collisions  of  the 
ions  and  neutrals.  Note  8  to  Table  1  explains  the 
mass-dependent  values  of  the  time  unit,  which  is  2.85  ps  for 
N2'*^,  and  Fp,  which  is  112  for  N2'^  when  njj  is  in  cm~^. 
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Figure  23.  Fraction  of  the  total  en  igy  in  the  ambient  ion  kinetic  energy 
after  461  time  units  as  a  function  of  the  number  density  of 
neutral  molecules.  The  lines  are  labeled  with  the  energy  of 
the  beam  ions  and  with  the  initial  fraction  of  ions  in  the 
beam  f{,Q.  The  solid  line  labeled  NO-IEL  represents  runs  where 
there  are  no  elastic  collisions  of  the  ions  and  neutrals.  Note 
8  to  Table  1  explains  the  mass-dependent  values  of  the  time 
unit,  which  is  2.85  ys  for  N2^,  and  Fp,  whicn  is  112  for  N2^ 
when  njyj  is  in  cm"^. 


For  the  case  where  £^5/!  is  two,  the  fraction  of  energy  going  into  the 
field  energy  decreased  from  0.029  to  0.011  as  the  number  density  increased 
from  about  3  x  10^*^  to  3  x  lO^^/cm^,  when  the  fraction  of  ions  that  are  in  the 
beam  is  0.25.  Presumably  this  decrease  in  the  fraction  of  the  energy  that 
goes  into  the  electric  field  is  the  result  of  elastic  collisions  of  the  ions 
and  the  electrons  with  the  neutral  molecules,  which  scatter  the  charged 
species  and  reduce  the  field  build-up.  However,  the  field  energy  is  rather 
noisy,  and  the  pressure  dependence  is  not  always  so  noticeable.  However,  the 
trend  is  for  the  fraction  of  the  energy  in  the  field  to  decrease  as  the 
pressure  increases. 

3. 5. 5. 4  Effects  of  Ion-Neutral  Elastic  Collisions.  Some  simulations  did  not 

include  the  effects  of  ion-neutral  elastic  collisions.  The  main  effects  of 
these  collisions  is  to  increase  the  thermal  energy  of  the  ambient  ions  and  to 
reduce  the  amount  of  CIV  ionization.  The  curves  labeled  NO-IEL  show  the 
results  when  no  ion-neutral  elastic  collisions  are  included.  Figure  23  shows 
the  greatest  change  as  the  fraction  of  energy  in  the  ambient  ions  is  much 
greater  when  the  ion  elastic  collisions  are  included.  The  increase  in  ambient 
ion  thermal  energy  is  to  be  expected,  as  we  are  colliding  fast  neutral 
molecules  with  slow  ambient  ions.  Eventually  the  increase  in  the  thermal 
velocity  of  the  ambient  ions  could  reduce  the  electron  heating,  as  the 

velocity  of  the  beam  ions  vjjj  >  Vgjj.^  criterion  for  the  modified 

two-stream  instability  (McBride,  et  al.,  1972).  However,  our  simulations  did 
not  heat  the  ambient  ions  this  much.  We  did  observe  that  the  electron 
ionization  was  some  30  percent  less  when  the  elastic  collisions  were  included, 
while  the  field  energy  was  nearly  the  same,  but  Table  1  shows  that  this  effect 
varies  with  nj^. 

3. 5. 5. 5  Molecular  Effects  of  Increased  Neutral  Density.  In  this  section  we 
discuss  the  effects  of  increased  nf^  on  the  number  of  electron  ionization 
collisions  and  electron  excitation  collisions,  as  well  as  the  effects  on  the 
various  types  of  elastic  collisions  and  charge  exchange  collisions.  These 
effects  can  be  understood  by  using  the  definitions  of  the  various  collision 
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frequencies  given  in  Eq.  (102).  Figures  24  through  28  show  plots  of  the 
results  for  the  first  1.2  ms.  Recall  that  there  initially  are  1024  electrons 
and  that  the  ions  are  distributed  initially  with  a  fraction  f^Q  in  the  beam. 
Thus  the  number  of  ion-neutral  elastic  collisions  of  ambient  ions,  plotted  in 
Fig.  24,  gives  two  curves,  one  for  each  f^Q.  Also  note  that  the  curvature  in 
this  plot  and  in  Fig.  27  results  from  the  depletion  of  the  ambient  ions.  In 
Fig.  27  there  are  fewer  charge  exchange  collisions  of  the  ambient  ions  after 
the  number  of  such  collisions  begins  to  approach  the  initial  number  of  ambient 
ions,  which  is  512  when  fjjQ  is  0.5  and  1019  when  fjjQ  is  0.005.  At  lower 
values  of  njq,  the  number  of  ambient  ion  collisions  in  Figs.  24  and  27  is 
approximately  twice  as  great  for  fjjo  =  0.005,  as  there  are  approximately  twice 
as  many  ambient  ions  as  for  fjjo  =  0.5.  In  Fig.  25,  the  elastic  collisions  of 
the  beam  ions  continue  to  grow;  the  slope  of  the  fjjo  =0.5  data  is  less  than 
that  for  the  f^o  =  0.005  data  because  the  relative  change  in  the  number  of 
beam  ions  is  less.  There  are  some  effects  of  changing  Eg^/I.  Figure  28  shows 
that  the  Eg^/I  =  3  data  are  better  at  slowing  the  ions  so  that  there  are  more 
resonant  charge  exchange  collisions,  and  even  a  few  more  ion-elastic 
collisions  as  shown  in  Fig.  25.  Figure  26  shows  that  the  number  of  electron 
elastic  collisions  does  depend  on  the  electron  heating.  Thus  at  larger  values 
of  fjjo  or  Egjj/I»  the  electrons  are  more  energetic,  so  that  a  larger  fraction 
of  them  have  energies  near  the  maximum  in  the  elastic  cross  section  near  2  eV. 

3. 5. 5. 6  Requirement  for  a  Minimum  Density  of  Beam  Ions.  The  discussion  of 
the  CIV  process  often  assumes  that  the  process  can  develop  from  a  few  initial 
seed  ions,  with  no  mention  of  any  requirement  of  a  minimum  number  of  seed 
ions.  Further  consideration  shows  that  CIV  ionization  will  not  occur  at  a 
reasonable  rate  unless  there  is  a  sufficient  number  density  of  beam  ions  so 
that  the  electrons  are  heated  enough  to  produce  an  reasonable  fraction  with  an 
energy  exceeding  the  ionization  potential.  That  is,  electron  heating 
can  occur  without  a  threshold,  but,  if  there  is  only  a  small  number  of  beam 
ions,  then  the  energy  supplied  to  the  electrons  by  the  slowing  of  the  fast 
ions  will  be  too  small  to  raise  the  energy  of  a  sufficient  number  of  electrons 
to  produce  any  ionization  events.  We  can  understand  this  minimum  fraction  of 
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ELASTIC  COLLISIONS  (1000) 


Figure  24.  Number  of  elastic  collisions  of  the  ambient  ions  with  neutrals 
during  461  time  units  as  a  function  of  the  number  density  of 
neutral  molecules.  The  lines  are  labeled  with  the  energy  of 
the  beam  ions  and  with  the  initial  fraction  of  ions  in  the 
beam  fjjo-  The  solid  line  labeled  NO-IEL  represents  runs  where 
there  are  no  elastic  collisions  of  the  ions  and  neutrals.  Note 
8  to  Table  1  explains  the  mass-dependent  values  of  the  time 
unit,  which  is  2.85  ys  for  N2''',  and  Fp,  which  is  112  for  N2^ 
when  n(^  is  in  cm“3.  The  initial  number  of  ambient  ions  is  1024 

(i-£bo)- 
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ELASTIC  COLLISIONS  (1000) 
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Figure  25.  Number  of  elastic  collisions  of  the  beam  ions  with  neutrals 
during  461  time  units  as  a  function  of  the  number  density  of 
neutral  molecules.  The  lines  are  labeled  with  the  energy  of 
the  beam  ions  Egj,  and  with  the  initial  fraction  of  ions  in  the 
beam  f^Q.  The  solid  line  labeled  NO-IEL  represents  runs  where 
there  are  no  elastic  collisions  of  the  ions  and  neutrals. 

Note  8  to  Table  1  explains  the  mass-dependent  values  of  the 
time  unit,  which  is  2.85  ps  for  N2^,  and  Fp,  which  is  112  for 
N2'*'  when  nj,}  is  in  cm~3.  xhe  initial  number  of  ambient  ions  is 
1024  (fbo). 
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Figure  26.  Number  of  elastic  collisions  of  the  electrons  with  neutrals  during 
461  time  units  as  a  function  of  the  number  density  of  neutral 
molecules.  The  lines  are  labeled  with  the  energy  of  the  beam  ions 
Eg|j  and  with  the  initial  fraction  of  ions  in  the  beam  f^o-  The 
solid  line  labeled  NO-IEL  represents  runs  where  there  are  no 
elastic  collisions  of  the  ions  and  neutrals.  Note  8  to  Table  1 
explains  the  mass-dependent  values  of  the  time  unit,  which  is  2.85 
ys  for  N2‘’',  and  Fp,  which  is  112  for  N2‘''  when  njq  is  in  cm~^.  The 
initial  number  of  electrons  is  1024. 
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Figure  27.  Number  of  charge  exchange  collisions  of  ambient  ions  with  neutrals 
during  461  time  units  as  a  function  of  the  number  density  of 
neutral  molecules.  The  lines  are  labeled  with  the  energy  of  the 
beam  ions  E|qjj  and  with  the  initial  fraction  of  ions  in  the  beam 
fjjQ.  The  solid  line  labeled  NO-IEL  represents  runs  where  there 
are  no  elastic  collisions  of  the  ions  and  neutrals.  Note  8  to 
Table  1  explains  the  mass-dependent  values  of  the  time  unit,  which 
is  2.85  ps  for  N2''’,  and  Fp,  which  is  112  for  N2'''  when  nj^  is  in 
cm"3.  The  initial  number  of  ambient  ions  is  1024  (l-fjjQ). 
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Figure  28.  Number  of  charge  exchange  collisions  of  beam  ions  with  neutrals 
during  461  time  units  as  a  function  of  the  number  density  of 
neutral  molecules.  The  lines  are  labeled  with  the  energy  of  the 
beam  ions  Egj,  and  with  the  initial  fraction  of  ions  in  the  beam 
fjjQ.  The  solid  line  labeled  NO-IEL  represents  runs  where  there 
are  no  elastic  collisions  of  the  ions  and  neutrals.  Note  8  to 
Table  1  explains  the  mass-dependent  values  of  the  time  unit, 
which  is  2.85  gs  for  N2^,  and  Fp,  which  is  112  for  N2'^  when  n]^ 
is  in  cm~3.  xhe  initial  number  of  beam  ions  is  1024  (fbo)- 
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the  ions  that  are  beam  ions,  by  considering  that  the  electron  heating  from  nj, 
fast  ions  that  initially  have  a  kinetic  “Energy  of  (=  mj,  will  be 

A  Eg  =  h  njj  Egjj  .  (105) 

Thus  the  average  increase  in  energy  is  h  nj,  Egb^^’^e^final*  Electrons  must  be 
energized  to  values  above  I  before  ionization  can  occur.  Thus  we  can  assume 
that  the  average  increase  in  energy  must  exceed  some  fraction  h  of  I,  and 
when  the  CIV  process  is  in  the  early  stages  with  (^’e^final  ~  then 


h  "b  Egb^'^e  >  hi  , 

(106) 

f  =  nb/ng  >  hl/(ri  Eg^j)  . 

(107) 

This  minimum  fraction  is  thus  a  function  of  h,  and  the  value  of  h  will  depend 
on  the  velocity  distribution  of  the  electrons.  If  this  distribution  is 
Maxwellian,  then  h  ~  0.20  is  a  reasonable  estimate.  The  electron  distribution 
is  expected  to  have  a  pronounced  tail  (Papadopoulus ,  1984).  When  the  drift 
velocity  of  the  beam  ions  is  small  (vj,  <  3  then  the  electrons  will  form 

a  tail  (Tanaka  and  Papadopoulos,  1983).  However  our  beam  ions  are  faster  than 
this  so  our  electron  distribution  may  not  have  as  pronounced  a  tail.  If  there 
is  a  significant  tail  in  our  simulations,  then  h  will  be  smaller. 

We  see  that  either  making  the  kinetic  energy  of  the  beam  greater  or  increasing 
the  efficiency  gives  the  condition  that  the  nj,  /  ng  ratio  can  be  smaller. 

In  our  simulations  with  no  charge  exchange  collisions  (and  also  with  no 
ion-neutral  elastic  collisions),  then  f,nin  ~  0.1  when  Eg|,  was  21  and  n  =  0.6. 
This  implies  that  h  =;  0.12  for  our  simulation  for  a  low  gas  density  around 
10l0/cm3. 

The  importance  of  this  effect  is  that  there  can  be  growth  of  a  plasma 
instability,  which  is  energizing  electrons,  but  the  tail  of  the  distribution 
must  have  electrons  with  energies  above  I  before  additional  ionization  can 
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occur.  Furthermore,  the  increase  of  electron  energy  will  be  at  least 
partially  distributed  among  the  bulk  of  the  electron  distribution,  which  leads 
to  a  requirement  of  f^in  CIV  ionization. 

3. 5. 5. 7  Effects  of  Kinetic  Energy  of  Fast  Beam  Ions.  When  the  kinetic  energy 
of  the  fast  beam  ions  is  increased,  the  electron  heating  increases.  As 
discussed  above,  the  increased  Eg^jj  of  the  fast  ions  allows  the  electrons  to  be 
heated  enough  to  produce  further  ionization  with  a  smaller  fraction  of  beam 
ions  than  would  be  required  if  the  Ej^jj  were  smaller.  Basically  the  increased 
energy  in  the  fast  ion  beam  produced  from  the  fast  neutrals  transfers  more 
energy  to  the  electrons,  which  makes  a  larger  fraction  of  the  electrons  have 
energies  above  I.  Figures  17  and  18  give  results  from  a  series  of  runs  where 
f|i}0>  the  fraction  of  ions  in  the  beam  initially,  was  either  0.005  or  0.5,  the 
ratio  of  the  Initial  kinetic  energy  of  the  beam  ions  to  the  ionization 
potential  varied  from  0.6  to  3,  and  the  number  density  varied  from  10^®  (q 
10l3  per  cm3.  Figures  19  and  20  show  the  rates  to  heat  0.001  of  the  electrons 
to  an  energy  of  1/2  or  I,  respectively.  Table  1  has  the  time  to  increase  the 
electron  density  by  a  fraction  of  0.001. 

3. 5. 5. 8  Importance  of  Charge  Exchange  as  an  Energy  Source  for  CIV.  Note  that 
Table  1  and  Figs.  17  through  20  show  that  it  is  possible  to  heat  electrons  to 
the  ionization  energy  even  when  the  kinetic  energy  of  the  beam  molecules  is 
less  than  the  ionization  energy.  Furthermore,  it  is  even  possible  to  produce 
an  increase  in  the  electron  density  when  this  is  true.  The  present 
simulations  have  not  followed  the  increase  over  a  long  period,  but  we  can 
anticipate  that  there  is  a  greater  increase  than  would  be  allowed  if  the  only 
energy  source  is  the  kinetic  energy  of  the  beam  ions  that  are  present 
initially  or  formed  from  ionization  by  fast  electrons.  The  reason  for  this 
greater  increase  is  that  the  resonant  charge  exchange  collisions  are  an 
additional  source  of  energy.  That  is,  if  each  beam  ion  transfers  about  3/5  of 
its  kinetic  energy  to  the  electrons  we  might  expect  that  there  would  be  an 
energy  threshold  of  (5/3)  I  (not  to  mention  the  fact  that  there  are  also 
losses  from  energy  transferred  to  excited  states  that  decay  before  they  were 
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ionized).  However,  there  are  two  sources  of  additional  energy:  charge 
exchange  with  ambient  ions  and  resonant  charge  exchange  of  slower  beam  ions 
with  faster  beam  molecules. 

Each  time  there  is  charge  exchange  with  an  ambient  ion,  there  is  an  increase 
in  the  kinetic  energy  in  the  ion  beam.  The  energy  increase  is  given  by 

^  ^cea  “  ^'''b  '^b^  ~  "'ia  '^ith^^'^^  »  (108) 

and  the  beam  velocity  is  faster  than  the  thermal  velocity  of  the  ambient  ions 
v^fh*  In  fact,  McBride  and  coworkers  (1972)  show  that  vjj  >>  v^^j^  is  the 
condition  for  the  modified  two-stream  instability  that  is  primarily 
responsible  for  heating  the  electrons  and  slowing  the  beam  ions.  Thus  the 
energy  available  to  heat  electrons  from  this  source  is  approximately  n^^  ri  Egj,, 
where  we  have  neglected  the  ambient  ion  thermal  energy.  Therefore  the  energy 
increase  per  electron  from  charge  exchange  with  ambient  ions  followed  by 
slowing  the  fast  ions  is  fg^e  *1  ^Kb»  where  f^^^g  is  the  fraction  of  ambient 
ions  that  have  charge-exchange  collisions  with  the  fast  neutrals.  This  can  be 
calculated  at  any  time  t  from  the  rate  constant  for  charge  exchange  kggg  as 

^ace  =  1  -  exp(-kcea  n^j  t)  .  (109) 

Thus  from  Eqs.  (107)  and  (109),  when  charge  exchange  with  ambient  ions  is  the 
dominant  means  of  forming  beam  ions,  CIV  ionization  will  occur  for 

Egb/I  >  h/(n  face)  •  (HO) 

If  we  also  consider  that  we  have  energy  losses  for  excitation,  then  we  expend 
an  average  energy  of  W  each  time  we  create  a  new  ion  pair,  where  W  ~  <I>  + 
K<X>.  Here  <I>  and  <X>  represent  the  average  ionization  and  excitation 
energy,  respectively.  If  we  let 

w  =  W  /  I  ,  (111) 
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we  can  take  j^'yantage  of  the  correlation  that  w  is  about  1.7  for  rare  gases 
and  2.1  tw  2.5  for  molecular  gases  for  fast  electrons  (Platzman,  1961).  Our 
tail  electrons  are  only  moderately  fast,  so  we  expect  our  w  values  to  be 
somewhat  higher.  However,  in  the  Platzman  analysis  there  is  a  contribution  of 
0.2  to  0.3  to  w  from  energy  that  goes  into  the  secondary  electrons,  and  this 
term  will  also  contribute  to  our  electron  heating.  Thus  we  may  expect  to  have 
w  values  around  1.5  to  3.  After  including  the  energy  loss  to  excitation,  the 
condition  for  CIV  is 

Egb/I  >  h  w/(n  face)  •  (112) 

For  example  if  h  =  0.12,  w  =  2.2,  and  t)  =  0.6,  then  ^  0.44/face-  Thus, 

Egb  !  I  values  as  low  as  0.44  could  have  some  CIV  ionization. 

In  order  to  increase  the  electron  density,  an  energy  of  A  ng  w  I  must  be 
transferred  to  the  tail  electrons.  However  we  get  an  increase  in  the  kinetic 
energy  of  the  beam  ions  from  each  new  ion  formed,  and  0  Ek^Jj  energy  will 
eventually  transferred  to  the  electrons.  Thus  the  energy  balance  gives  the 
condition  for  CIV  to  grow  as 

Egb/I  >  w  (h+5)/(n(face  +  S) 1  *  (113) 

where  6  =  A  ng/(ng)Q  is  the  ratio  of  the  increase  in  electron  density  to  the 
initial  electron  density.  Using  the  same  values  in  this  equation  gives 
Egij/I  >  3.7  (0.12  +  5)/(face  +  *)«  2. 3/(0. 5  +  fggg)  for  6  =  0.5.  Thus 

Eg^b^l  ^  1-5  for  this  example.  Even  if  h  is  as  low  as  0.01,  w  is  as  low  as 
1.5,  t)  is  as  large  as  0.8,  then  Eg^/I  >  0.96/(face  0.5)  when  &  =  0.5.  Thus 

Egb^I  ^  1  even  for  this  example  when  face  ^  0.46.  Thus  it  is  possible  to  meet 
the  Alfven  criterion  of  =  1  for  this  example.  Moreover  we  still  have  an 

additional  energy  source  available  when  resonant  charge  exchange  is  rapid. 

The  energy  gain  from  each  resonant  charge  exchange  is 

(^  f-Kblcer  ~  '’’b('^b^  “  '^bs^)22  ,  (114) 
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where  vjjg  is  the  velocity  of  the  beam  ion  that  has  been  slowed.  If  we  take 
Vbs  =  Vb/2,  then 


EKb)cer  =  0. 75(mbVb2/2)  .  (115) 

Note  that  there  can  be  many  resonant  charge  transfer  collisions.  Thus  the 
total  kinetic  energy  in  the  directed  beam  can  actually  increase  if  the 
resonant  charge  exchange  frequency  “LH  >  that  the  resonant  charge 

exchange  feeds  energy  into  the  beam  faster  than  the  plasma  instabilities  can 
remove  it.  The  equilibrium  condition  would  be  when  the  charge  exchange  rate 
is  equal  to  the  ion-slowing  rate,  which  is  the  electron  heating  rate.  This 
rate  will  be  related  to  the  modified  two-stream  instability  growth  rate  y, 
which  is  estimated  as 

Y  ~  0.6  fO-^  (OLH  .  (116) 

where  <»)lh  1^  ^he  lower  hybrid  frequency  (Akimoto,  et  al.,  1985).  If  we  assume 
that  the  resonant  charge  exchange  frequency  is  equal  to  y,  then 

^cer  =  <<^cer  Vr>  :  0.6  fO-^  «lh  » 

where  the  relative  velocity  of  the  ion-neutral  collision  Vj.  ~  vjj  /  2.  For 
example,  for  N2,  <Jcer  ~  ^  lO'l^  (Stebbings,  et  al.,  1963),  so  for 

Vj.  =  4  X  10^  cm/s  and  wlh  =  Wj/fl  +  =  2  x  10^  s“l  when  the  ratio 

of  the  electron  plasma  frequency  to  the  electron  cyclotron  frequency  Qg  is 
2.26,  which  is  the  value  corresponding  to  =  10^/cm^  and  a  magnetic  field  of 
0.45  Gauss,  njj  =  0.9  x  lol^  f-0.4  _  (0.9  to  2)  x  lol^  cm  3  for  f  values  from 
0.1  to  1.  Thus  at  gas  densities  above  this  it  should  be  possible  to  maintain 
the  full  energy  of  the  fast  ion  beam  through  resonant  charge  transfer. 

Figure  22  shows  that  our  simulations  meet  this  density  condition  for  nj^  around 
I0I2  cm~3.  Actually  there  will  also  be  elastic  ion-neutral  collisions,  and 
these  will  change  the  requirement.  However  there  should  be  many  gases  that 
have  CTcgj-  values  large  enough  that  the  kinetic  energy  of  the  beam  ion  is 
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replenisned  by  resonant  charge  exchange  faster  than  it  is  transferred  to  heat 
the  plasma  electrons. 

In  one  series  of  runs  we  compared  the  effects  of  including  resonant  charge 
exchange  and  elastic  collisions  of  the  beam  ions  with  the  neutral  beam 
molecules.  We  compared  three  runs:  with  both  effects  (run  86),  with  neither 
effect  (88),  and  with  only  the  ion-neutral  collisions  (91).  At  a  neutral 
density  around  lO^^/cm^  and  Egj,  =  I,  we  find  that  the  resonant  charge  exchange 
enhancement  is  less  than  the  ion-neutral  reduction.  Thus  there  is  more  CIV 
ionization  for  the  case  with  neither  effect.  This  seems  to  result  from  the 
greater  build-up  of  the  electric  field  energy  when  there  are  no  ion-neutral 
collisions.  As  a  result  of  the  lower  field,  the  kinetic  energy  of  the  fast 
beam  builds  up  more  and  a  smaller  fraction  of  this  energy  is  transferred  to 
electron  kinetic  energy.  When  both  processes  are  present  there  is  an  increase 
in  the  drift  energy  of  the  beam  ions  by  over  40  percent,  but  the  electron 
kinetic  energy  is  almost  30  percent  smaller. 

3.5.6  Summary  of  Simulation  Results 

We  have  shown  that  for  conditions  present  in  the  high-density  region  of  a  beam 
of  fast  gas  molecules  that  charge  exchange  and  elastic  collisions  are 
important  physical  phenomena.  It  is  likely  that  the  experimentally  observed 
onset  of  CIV  effects  at  v^.  calculated  from  the  ionization  energy  is  the  result 
of  charge  exchange  effects  that  feed  additional  energy  into  the  process  to 
compensate  for  the  losses  from  the  failure  of  the  fast  ions  to  feed 
100  percent  of  their  energy  to  the  plasma  electrons  (i.e.,  h  <  1)  and  from  the 
losses  to  form  excited  states  rather  than  ions  (w  >  1).  When  the  fraction  of 
beam  ions  is  low,  the  electrons  will  heat,  but  not  enough  to  produce  further 
ionization.  Charge  exchange  with  ambient  (background)  ions  will  produce  the 
initial  fast  ions  when  the  beam  number  density  is  large  enough  and  the  beam 
molecules  have  an  ionization  potential  low  enough  that  this  process  is 
exothermic.  Some  charge  exchange  has  to  take  place  before  the  electrons  will 
heat  effectively.  This  imposes  low  pressure  limits  on  the  gas  beam.  At 
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higher  beam  gas  densities,  the  drift  energy  of  the  beam  ions  will  increase 
until  the  electric  field  energy  is  maintained  high  enough  to  heat  electrons  in 
spite  of  the  disrupting  effects  of  elastic  collisions.  Resonant  charge 
exchange  will  be  less  important  in  increasing  the  energy  of  the  beam  ions  when 
dimer  ions  are  present,  unless  the  dimer  concentration  is  a  large  fraction  of 
the  molecules  in  the  beam.  This  may  reduce  the  importance  of  dimer  (and 
higher  n-mer)  formation  in  the  CIV  process.  However,  the  effect  of  resonant 
charge  exchange  is  to  lower  the  energy  requirement  for  the  fast  ion  beam,  so 
any  contribution  to  the  beam  from  dimers  will  be  an  additional  way  to  produce 
CIV  ionization  with  a  low  energy  beam. 
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APPENDIX 


This  appendix  includes  the  as  delivered  Phase  0/1  Accident  Risk  Assessment 
Report  for  the  proposed  experiment  to  study  the  role  of  the  critical  ioniza¬ 
tion  velocity  effect  in  producing  enhanced  plasmas  around  the  Space  Shuttle. 
The  experiment  design  has  continued  to  evolve  since  this  doucment  was  deliv¬ 
ered  and  some  of  the  information  provided  within  is  no  longer  appropriate. 

For  example,  the  test  gases  are  now  specified  as  neon,  xenon,  nitric  oxide  and 
carbon  dioxide,  and  the  remote  plasma  probe  is  now  re-designed  as  a  passive 
receiver.  The  present  experiment  concept  is  described  in  Volume  II  of  this 
report,  which  includes  all  viewgraphs  presented  at  the  CIV  experiment 
Preliminary  Design  Review. 

The  safety  issues  discussed  below  remain  appropriate  in  any  event. 
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1 .  INTRODUCTION 


This  report  provides  an  Accident  Risk  Assessment  Report  (ARAR) ,  equiv¬ 
alent  to  a  Phase  I  Safety  Data  Package,  for  the  Critical  Ionization  Velocity 
(CIV)  Experiment.  It  identifies  and  evaluates  potential  hazards  for  equipment 
and  operations  required  for  flight  operations  of  the  CIV  equipment. 

This  report,  in  compliance  with  JSC  13830A,  provides  information  neces¬ 
sary  to  verify  that  design  procedures  and  plans  comply  with  NHB  1700. 7A.  It 
should  be  noted  that  some  data  concerning  the  payload  are  still  to  be  deter¬ 
mined  (TBD)  and  this  document  will  be  updated  as  information  is  available. 

The  CIV  experiment  is  a  part  of  the  Infrared  Background  Signature  Survey 
(IBSS)  mission  sponsored  by  the  Strategic  Defense  Initiative  Organization 
(SDIO).  The  CIV  experiment  will  be  designed,  developed,  and  fabricated  by 
Physical  Sciences  Inc.  (PSI)  under  contract  to  the  Air  Force  Geophysical 
Laboratory  (AFGL).  AFGL  is  developing  and  providing  the  CIV  experiment  with 
the  Payload  Support  Subsystem  as  described  in  Subsection  4.4. 

The  CIV  experiment  involves  the  sequential  release  of  four  gases  from  the 
Orbiter  payload  bay  and  the  observation  of  the  resulting  gas  cloud  from  the 
IBSS  system  (IBSS,  AIS,  and  TV  cameras  mounted  on  the  RMS  deployed  or  free 
flying  SPAS  platform)  and  from  a  CIV  monitor  subsystem  in  the  payload  bay. 

The  CIV  experiment  consists  of  three  subsystems:  1)  a  gas  release  subsystem 
with  four  pressure  vessels  and  pneumatic  control  hardware;  2)  a  monitor  sub¬ 
system  with  a  series  of  radiometers,  a  Langmuir  capacitance  probe,  and  a 
remote  plasma  monitor  utilizing  a  low  power  radio  frequency  transmitter/ 
receiver;  and  3)  a  payload  support  subsystem  with  a  payload  controller,  data 
acquisition  electronics,  data  tape  recorder,  and  a  battery. 

The  CIV  experiment  is  mechanically  configured  to  mount  on  two  of 
Rockwell's  Adaptive  Payload  Carriers  (APC)  attached  to  the  Orbiter  payload  bay 
longeron.  The  electrical  interface  with  the  orbiter  is  through  five  switches 
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of  the  Payload  Standard  Switch  Panel  (PSSP)  assigned  to  the  IBSS  mission  which 
provide  power  to  activate  relays  and  solenoid  valves  that  control  CIV  func¬ 
tions.  The  major  power  source  for  CIV  equipment  is  a  self-contained  battery 
pack. 
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2.  SAFETY  CONSIDERATIONS 


2 . 1  Safety  Certification 


AFGL  and  SDIO  management  approval  and  submittal  of  this  ARAR  for  the 
safety  review  process  constitutes  certification  of  completeness  in  accordance 
with  the  requirement  of  JSC  13830  for  Phase  0/1  Safety  Data  Package. 

2.2  Program  Safety  Status  Summary 

Development  of  this  ARAR  and  hazard  analysis  is  at  the  Phase  I  stage. 

CIV  system  description,  drawings,  schematics,  and  operation  have  been 
sufficiently  defined  to  allow  assessment  of  the  hazards  relating  to  the 
payload. 

2.3  Non-Compliant  Items 

Neither  waivers  nor  deviations  of  safety  requirements  for  CIV  are  antici¬ 
pated  at  this  time. 

2 .4  Safety  Relevant  Accident/Incidents 

None  at  this  time. 

2.5  Hazardous  Materials 


No  hazardous  materials  are  used  in  CIV  except  for  a  lithium-iodine  power 
cell  utilized  in  AFGL's  Payload  Support  System  (refer  to  attachment). 
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SYSTEM  DESCRIPTION 


3  . 


The  Critical  Ionization  Velocity  (CIV)  Experiment  will  investigate  a 
physical  mechanism  for  creating  ionized  species  whereby  neutral  molecules 
injected  into  the  ambient  low-earth  atmosphere  at  sufficient  velocity  and  in 
the  pressure  of  the  earth’s  magnetic  field  may  ionize.  The  CIV  experiment 
will  test  this  hypothesis  by  releasing  four  different  gases  from  the  orbiter 
payload  bay  and  monitoring  the  gas  cloud  from  a  self-contained  CIV  diagnostic 
package  and  from  the  free  flying  (or  RMS  deployed)  IBSS  platform. 

3 . 1  Critical  Ionization  Velocity  Experiment 

The  CIV  experiment  shown  in  Figure  1  consists  of  three  subsystems;  the 
gas  release  subsystem,  the  monitor  subsystem,  and  the  payload  support 
subsystem. 

The  gas  release  subsystem  consists  of  a  pressure  vessel,  safety  relief 
valve,  a  pressure  regulator,  dual  solenoid  valves  and  exit  nozzle  for  each  of 
four  test  gases.  The  expansion  nozzle  will  direct  released  gases  parallel  to 
the  orbiter  +Z  axis.  Pressure  and  temperature  sensors  monitor  the  gas  condi¬ 
tions  in  each  pressure  vessel.  The  four  test  gases  are  argon  (Ar) ,  carbon 
dioxide  (CO2) ,  nitrous  oxide  (NO),  and  TBD. 

The  monitor  subsystem  contains  three  diagnostic  instruments:  a  series  of 
bandpass  filtered  radiometers,  a  Langmuir  capacitance  probe  and  a  remote 
plasma  monitor.  The  radiometers  are  calibrated  silicon  or  PbS  photodiodes 
with  bandpass  interference  filters  and  apertures  to  limit  the  field  of  view. 
The  Langmuir  probe  is  a  small  conducting  sphere  mounted  on  a  support  rod  which 
is  used  to  make  a  passive  measurement  of  the  capacitance  of  the  space  around 
the  probe.  The  remote  plasma  monitor  uses  a  10  watt  radio  frequency  trans¬ 
mitter  and  receiver  to  measure  the  reflective  properties  of  the  ionized  gas 
cloud  as  the  radio  frequency  is  varied. 
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Figure  1.  CIV  Experiment  Configuration 


The  payload  support  subsystem  contains  a  programmable  payload  controller 
that  coordinates  CIV  experiment  operations,  including  activation  of  the 
solenoid  valves  for  gas  release  in  a  preprogrammed  sequence  (additional  PSSP 
switches  can  optionally  activate  solenoid  valves  directly  to  override  the 
preprogrammed  sequence).  The  data  acquisition  module  multiplexes  and  formats 
data  from  the  monitor  subsystem  instruraents  and  the  gas  release  diagnostic 
sensors  and  directs  the  data  to  the  tape  recorder.  The  payload  support 
subsystem  also  contains  a  28-volt,  50-amp-hr  silver  zinc  battery  and  DC/DC 
converters  to  supply  +5  and  ±5  volts  to  payload  electronics. 

The  gas  release  subsystem  is  mounted  on  one  APC  and  the  monitor  subsystem 
and  payload  support  subsystem  are  mounted  on  the  second  APC.  A  wire  harness 
will  provide  electrical  connections  for  power  and  control  between  the  two  APC 
which  may  be  mounted  on  opposite  sides  of  the  payload  bay. 

3 .2  Flight  Operations 

Operation  of  CIV  functions  are  initiated  by  a  crew  member  on  the  aft 
flight  deck  through  the  PSSP.  Three  switches  are  utilized  for  baseline  CIV 
operations,  and  two  additional  switches  allow  optional  functions  to  be  accom¬ 
plished.  The  switch  schematic  diagrams  are  shown  in  Figure  2.  The  particular 
switch  numbers  referred  to  are  for  information  only.  Any  equivalent  switch 
would  be  acceptable. 

CIV  remains  in  its  unpowered  state  until  shortly  before  a  planned  gas 
release  when  the  crew  member  activates  switch  S7  which  drives  a  relay  that 
applies  CIV  battery  power  to  the  payload  controller.  The  payload  controller 
immediately  applies  power  to  the  radiometers,  the  Langmuir  probe,  and  the  data 
recorder.  The  payload  controller  also  directs  power  to  the  remote  plasma 
monitor  transmitter;  however,  dual  relays  interrupt  this  power  to  inhibit 
transmitter  operation  unless  switch  S3  is  activated.  The  activation  of  switch 
S3  by  the  crew  member  completes  the  power  up  and  initialization  of  the  CIV. 
With  S7  and  S3  both  activated,  the  monitor  instruments  are  operating  and  data 
is  continually  recorded. 
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Figure  2.  CIV  Control  Schematic 


Prior  to  gas  release  appropriate  orbiter  attitude,  IBBS  attitude,  and 
IBSS  operational  mode  must  be  achieved.  Then  gas  releases  are  initiated  in 
one  of  two  methods.  For  baseline  gas  releases,  switch  S8  is  activated  which 
produces  an  interrupt  on  the  payload  controller  and  a  preprogrammed  sequence 
of  events  begins.  In  the  sequence  the  solenoid  valve  for  each  test  gas  is 
opened  in  turn  for  a  10-second  gas  release  with  a  null  period  between  each  gas 
release  of  approximately  30  seconds. 

The  optional  gas  release  method  uses  switches  S1 1  and  S12  to  operate  each 
solenoid  valve  directly,  bypassing  the  payload  controller  sequencing.  For 
such  manually  controlled  gas  releases,  the  crew  member  would  first  position 
switch  S1 1  to  the  proper  position,  then  when  the  gas  release  is  desired 
operate  switch  S12.  Table  1  gives  typical  switch  positions  for  various  gases. 

Table  1 .  Optional  Gas  Release  Control  Switches 


S11  Position 

SI 2  Position 

Test  Gas  Released 

Up 

Up 

A 

Up 

Down 

B 

Down 

Up 

C 

Down 

Down 

D 

For  the  optional  control  method,  the  crew  member  would  be  required  to  manually 
time  each  gas  release  and  the  null  time  between  gas  releases. 

After  the  CIV  measurement  period  is  over,  switches  S3  and  S7  are  deacti¬ 
vated  to  terminate  data  acquisition  and  remove  power  from  the  CIV  payload. 


Five  sets  of  gas  release  sequences  are  currently  planned  with  the  orbiter 
and  IBSS  in  various  attitudes  and  configurations.  Depending  on  the  length  of 
time  between  gas  releases,  the  crew  member  could  remove  power  from  CIV  or 
simply  leave  the  monitor  subsystem  operating. 

A  preliminary  timeline  for  proposed  CIV  operations  is  presented  in 
Figure  3. 
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Figure  3.  CIV  Flight  Operations 
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4.  SYSTEM  SAFETY  CONSIDERATIONS 


The  CIV  experiment  system  has  been  analyzed  to  identify  safety  critical 
components  and  functions  in  accordance  with  NHB  1700. 7A.  This  section 
describes  the  results  of  hazard  identification  and  discusses  the  general 
approach  to  hazard  control  in  the  CIV  subsystems. 

4 . 1  Payload  Safety  Requirements  Applicability  Matrix 

The  hazard  matrix  for  the  CIV  experiment  is  presented  in  Table  2. 

4 .2  Gas  Release  Subsystem 


Preliminary  characteristics  of  the  gas  release  subsystem  are  presented  in 
Table  3.  Identified  hazards  include  pressure  vessel  explosion  and  collision 
of  failed  structural  components. 

The  pressure  vessels  will  be  purchased  from  a  manufacturer  who  has  pre¬ 
viously  qualified  the  pressure  vessel  design  for  space.  As  such,  pressure 
vessels  will  be  designed,  built,  and  tested  in  accordance  with  MIL-STD- 1522A 
or  NSS/HP  1740.1.  The  specified  safety  factor  is  1.5  with  a  leak-before-burst 
design.  Redundant  pressure  relief  devices  (safety  relief  valve  and  burst 
disk)  are  incorporated  into  the  gas  control  system  to  prevent  overpressuriza¬ 
tion  under  all  thermal  conditions. 

All  pressurized  components  including  lines,  fittings,  solenoid  valves, 
and  regulators  will  be  space  qualified  or  military  specification  parts  to  meet 
guidelines  of  MIL-STD- 1 522A  or  NSS/HP  1740.01.  All  pressurizea  components 
will  be  rated  for  maximum  operating  pressure. 

The  four  test  gases  are  argon  (Ar) ,  carbon  dioxide  (CO2) ,  nitrous  oxide 
(NO),  and  TBD.  All  are  non-hazardous  for  all  flight  operations.  The  thrust 
produced  by  the  CIV  gas  release  is  small  compared  to  the  vernier  reaction 
control  system  (VRCS)  thrusters.  Approximately  5N  (1.2  Ibp)  of  thrust  is 


124 


Table  3.  Gas  Release  Subsystem  Operating  Characteristics 


Cylindrical  pressure  vessels 

Volume 

- 

14  liter 

Length 

- 

60  cm 

OD 

- 

20  cm 

Weight 

- 

1 4  kg 

Operating  pressure  - 

125  bar 

CO2 

NO 

Ar 

TBD 

Total  mass  of  gas  (kg) 

3.0 

2.3 

2.8 

TBD 

Fill  pressure  (bar) 

35.0 

125.0 

4  125.0 

Gas  release  rate  (g/s) 

11.0 

7.5 

TBD 

Typical  release  time  (s) 

5.0 

5.0 

B 

5.0 

produced  by  the  5-second  gas  release  and  roughly  50  grams  of  gas  is  released. 
By  comparison,  the  vernier  reaction  control  system  (VRCS)  thrusters  produced 
1 1 4N  (25  Ibp)  of  thrust  and  produce  a  higher  mass  flux  of  species  including 
H2O,  ^2 *  ^^2*  CO,  and  CO2. 

4.3  Monitor  Subsystem 

The  monitor  subsystem  consists  of  three  Instruments  for  diagnostic  moni¬ 
toring  of  the  gas  cloud:  photodiode  arrays,  a  Langmuir  capacitance  probe,  and 
a  radio  frequency  transmitter/receiver.  The  photodiode  arrays  are  passive 
filter /detector  packages  which  record  the  radiance  of  the  gas  cloud  in  various 
wavelength  bands.  The  Langmuir  capacitance  probe  is  a  6-cm  steel  sphere  which 
is  extended  on  a  support  rod  approximately  1  meter  above  the  APC  top  surface. 
The  remote  plasma  monitor  uses  a  10-watt  transmitter  to  sweep  the  radio  fre¬ 
quency  while  monitoring  the  reflected  signal  with  a  receiver.  The  hazards 
identified  for  this  subsystem  are:  radiation,  and  collision  of  failed  struc¬ 
tural  elements. 
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The  active  RF  plasma  monitor  will  be  designed  to  meet  the  EMI  allowable 
field  strength  limits  specified  in  ICD  19001  and  an  analysis  will  be  performed 
to  insure  that  the  RF  radiation  is  completely  compatible  with  all  payload  bay 
door  open  operations  of  the  orbiter  and  other  payloads.  The  transmitting 
antenna  will  be  directed  in  the  +Z  direction  through  the  payload  bay  doors. 

The  plasma  monitor  will  only  be  operated  during  each  of  four  to  six  test  gas 
release  sequences  and  several  short  periods  prior  to  the  releases  to  measure 
background  levels.  The  subsystem  power  will  be  off  during  launch  and  landing. 
Two  additional  relays  are  used  as  inhibits  to  prevent  inadvertent  o^jeration  of 
the  transmitter  due  to  payload  controller  EMI  susceptability  or  a  relay  fail¬ 
ure.  Because  the  IBSS  system  is  expected  to  be  free  flying  during  plasma 
monitor  operation,  the  appropriate  IBSS  specification/requirement  (TBD)  for 
allowable  radiated  field  strength  will  be  observed. 

4 .4  Payload  Support  Subsystem 

The  payload  support  subsystem  is  an  Air  Force  Geophysics  Laboratory 
(AFGL)  package  that  provides  instrument  control,  data  recording,  and  battery 
power.  This  subsystem  was  conceived  to  be  general  purpose  flight  hardware  for 
small  scientific  payloads.  AFGL  is  currently  developing  the  subsystem  for  the 
Visual  Photometric  Experiment  (VIPER)  which  is  mounted  in  a  GAS  cannister  with 
a  motorized  door  assembly.  The  VIPER  program  has  completed  the  Phase  I  safety 
review  and  the  data  package  is  attached  as  Appendix  A.  Currently  the  payload 
support  subsystem  for  the  VIPER  program  is  built  and  component  tested  with 
subsystem  level  tests  ongoing.  Hazards  associated  with  this  subsystem  are: 
explosion,  corrosion,  and/or  contamination  due  to  battery  failure,  fire,  and 
collision  of  failed  structures. 

Explosion  and  corrosion  are  hazards  associated  with  the  main  silver-zinc 
battery  used  to  provide  primary  power.  These  hazards  will  be  controlled  by 
housing  this  battery  in  a  sealed,  corrosion  proof,  nylon/fiberglass  box.  Eagle 
Picher  (EP)  Part  No.  16A10429-1.  This  is  vented  to  the  NASA  supplied  pressure 
relief  valves.  This  case  will  be  pressure  tested  to  22  psi.  The  free  volume 
of  the  box  will  be  minimized  to  reduce  excess  accumulation  of  gas  and  No.  22 
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Fiberglass  B-50,  aianufacturecl  by  Johns  Mansville,  will  be  installed  to  contain 
any  escaping  electrolyte. 

A  second  battery  is  located  within  the  payload  controller  to  power  a 
real-time  clock  when  the  main  battery  is  off.  This  is  a  Catalyst  Research 
Corp.  Model  3440  solid  state  lithium  iodine  2.8  volt  cell  with  a  capacity  of 
650  milliamp-hours.  Its  total  weight  is  17.8  grams  with  lithium  weight  of 
0.247  gram.  It  is  soldered  directly  to  the  controller  PC  board  and  its 
hermetically  sealed,  welded  nickel  case  needs  no  venting.  Because  of  the 
battery's  small  size,  no  hazard  is  associated  with  it. 

Because  the  payload  support  subsystem  will  be  mounted  directly  to  an  APC 
and  not  in  a  GAS  cannister,  the  electronics  housings  will  be  modified  from 
their  current  design  to  make  sealed  enclosures  that  fully  satisfy  the  appro¬ 
priate  structural  and  containment  requirements. 

4.5  Structural  Mount  Subsystem 

The  structural  mount  subsystem  consists  of  two  APC  plates  and  mounting 
hardware  for  each  payload  component.  Each  electrical  component  will  be 
installed  in  its  own  structural  housing  and  bolted  directly  to  the  APC.  A 

mounting  bracket  and  support  structure  will  be  designed  to  attach  the  four 

pressure  vessels  and  associated  fittings  to  the  APC. 

The  structural  design  will  comply  with  NHB  1700. 7A,  Para.  802.  Newly 

designed  structural  components  will  utilize  a  factor  of  safety  greater  than  or 
equal  to  1.4.  Fracture  analysis  will  be  performed  on  critical  elements  and 
redundant  or  fail-safe  designs  will  be  incorporated  where  possible.  Design 
load  factors  will  be  based  on  Information  from  ICD  19001  and  from  previous 
analyses  for  qualification  of  the  APC  structures. 

The  control  and  test  plan  for  structural  verification  is  presently  being 
deve loped. 


5 .  SAFETY  DATA 


This  section  incluJes  hazard  reports  for  each  identified  hazard 
associated  with  the  CIV  experiment  flight  operations. 


5 . 1  Hazard  List  (Flight  Operations) 


Equipment/ 

Subsystem 

Hazard 

Report 

# 

Title 

Applicable 

Paragraph 

NHB  1700.7 

Status 

All 

FI 

Collision  of  failed 
structural  elements 

208 

Open 

Gas  release 
subsystem 

F2 

Failure  of  pressurized 
components 

208 

Open 

Electrical 

F3 

Ignition  of  flammable 
atmospheres  or  materials 

209,213,219 

Open 

Electrical 

F4 

Exposure  of  STS  to  EMI  or 
transmitter  radiation 

212.2 

Open 

Electrical 

F5 

Battery  explosion  or 
release  of  corrosive 
materials  into  payload  bay 

209 

Open 

5.2  Hazard  Reports 


The  hazard  reports  and  appropriate  annexes  follow. 
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PAYLOAD  HAZARD  REPORT 


NUMBER:  FI 


PAYLOAD:  CIV 

PHASE:  0/1 

SUBSYSTEM:  Structures 

Collision 
HAZARD  GROUP:  Corrosion 

DATE :  Aug  87 

HAZARD  TITLE: 

Collision 

of  Failed  Structural  Elements 

APPLICABLE  SAFETY  REQUIREMENTS:  NHB 

1700.7,  Para.  208 

HAZARD  CATEGORY: 

X  Catastrophic 
Critical 

DESCRIPTION  OF 

HAZARD : 

Collision  of 

failed  CIV 

parts  with 

the  orbiter  and/or 

other  payloads 

HAZARD  CAUSES: 


Structural  failure  of  CIV  components  due  to: 

1.  Inadequate  structural  design 

2.  Faulty  manufacturing  and/or  fabrication  techniques 

3.  Stress  corrosion  cracking 


HAZARD  CONTROLS: 

J.  Use  IBSS  structural  design  criteria: 

a.  Ultimate  factor  of  safety  >  1.4  to  IBSS  design  limit  loads 

b.  Fail  safe  design  (redundant  structures)  where  possible 

c.  Fracture  control  to  safe  life  >  4  times  mission  life  where  redundant 

structures  not  possible 

2.  Use  safe  design  practice  as  in  1. 

3.  Ail  materials  to  be  selected  in  accordance  with  MSFC-SPEC-522  and  all 
cleaning  agents  to  be  controlled  by  specification  and  manufacturing 
logs . 


SAFETY  VERIFICATION  METHODS: 

Combination  of  analysis,  inspection,  and  test 


STATUS : 

Open  _ 

Approval _ Payload  Organization _ STS 

Phase  1 _ _ 


Phase  2 


Phase  3 
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PAYLOAD  HAZARD 

REPORT 

NUMBER:  F2 

PAYLOAD:  CIV 

PHASE ;  0/ 1 

SUBSYSTEM:  Pressure  Systems 

HAZARD 

GROUP ; 

Explosion 

DATE:  Aug  87 

HAZARD  TITLE: 

Failure  of  Pressurized 

Components 

APPLICABLE  SAFETY  REQUIREMENTS: 

NHB 

1700. 7A, 

Para.  208 

HAZARD  CATEGORY: 

X  Catastrophic 
Critical 

DESCRIPTION  OF 

HAZARD: 

Explosion  of  pressure  vessels,  tubing  and/or  pressurized  components 
causing  release  of  fragments  into  the  payload  bay 


HAZARD  CAUSES : 

1 .  Inadequate  design 

2.  Overpressurization  due  to  excessive  temperature  environment 


HAZARD  CONTROLS: 

1.  Pressure  vessels  designed,  manufactured,  and  tested  in  accordance 
with  MIL-STD-1522  or  NSS/HP  1740.1  with  safety  of  1.5  or  greater. 

All  other  pressurized  components  will  be  space  qualified  and/or 
military  specification  hardware  with  ultimate  factor  of  safety  _>  4.0. 

2.  Space  qualified  pressure  relief  mechanism  will  vent  gas  in  case 
of  overpressurization  under  all  thermal  conditions.  See  attached 
schematic  of  pneumatic  system  layout. 


SAFETY  VERIFICATION  METHODS: 

1 .  Analysis  and  environmental  test  under  full  pressure  -  pressure 
vessel  log 

2.  Test  to  verify  safety  relief  on  overpressure 
STATUS: 

OPg-» _ 

Approval _ Payload  Organization _ STS 

Phase  1 _ 

Phase  2 _ 

Phase  3  _  _ 
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SAFETY 

RELIEF 

VALVE 


CIV  Gas  Release  Subsystem  Schematic 
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PAYLOAD  HAZARD  REPORT 


NUMBER:  F3 


PAYLOAD: 

CIV 

PHASE: 

0/1 

SUBSYSTEM: 

All 

HAZARD  GROUP:  Explosion,  Fire 

DATE: 

Aug  87 

HAZARD  TITLE:  Ignition  of  Flammable  Atmospheres  or  Materials 


APPLICABLE  SAFETY  REQUIREMENTS:  NHB  1700.7A,  HAZARD  CATEGORY: 

Para.  208,213,219  X  Catastrophic 

Cri tical 


DESCRIPTION  OF  HAZARD: 

Failure  and/or  overheating  of  electrical  components  causes  ignition  of 
flammable  atmospheres  or  materials 


HAZARD  CAUSES: 

1 ,  Electrical  ignition  sources 

2.  Unrestricted  flame  propagation  paths 


HAZARD  CONTROLS: 

1.  Proper  design  of  electrical  circuitry,  components,  and  cables.  Fusing 
utilized  to  limit  current  to  appropriate  level  for  wire  size. 

2.  Exposed  material  selection  in  accordance  with  NHB  8060. IB  to  minimize 
flame  propagation.  See  attached  wiring  diagrams. 


SAFETY  VERIFICATION  METHODS: 

-  Analysis 

-  Review  of  material  lists  and  materials  assessment 


STATUS : 
Open 


Approval 
Phase  1 


Phase  2 


Phase  3 


Payload  Organization 


STS 
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A-6295 


Payload  Support  Subsystem  Power  Distribution  Schematic 


PAYLOAD  HAZARD  REPORT 

NUMBER:  F4 

PAYLOAD:  CIV 

PHASE 

0/1 

SUBSYSTEM:  Electrical  HAZARD  GROUP:  Radiation 

DATE: 

Aug  37 

HAZARD  TITLE:  Exposure  of  STS  to  EMI  or  Transponder 

Radiation 

APPLICABLE  SAFETY  REQUIREMENTS:  NHB  1700.7A, 

Para.  212.2 

HAZARD 

X 

CATEGORY : 
Catastrophic 
Critical 

DESCRIPTION  OF  HAZARD: 


1 .  Excessive  electromagnetic  radiation  of  CIV  electronics  interferes  with 
orbiter  or  other  payloads 

2.  Radio  wave  monitor  transponder  generates  excessive  interference  with 
orbiter  or  other  payloads 


HAZARD  CAUSES: 

1.  EMI  levels  exceed  allowance  of  JSC  07700  VOL  XIV  ATT  1  due  to  improper 
design  and/or  faulty  shielding 

2a.  Radiation  level  from  radio  wave  monitor  exceeds  level  specified 
2b.  Inadvertent  turn  on  of  transmitter  when  bay  doors  are  closed. 


HAZARD  CONTROLS: 

la.  Design  to  meet  ICD  19001  requirements 

lb.  Electrical  bonding  compliant  with  MIL-B-5078B/ICD  19001 

2a.  Radio  wave  monitor  designed  to  meet  JSC  07700  VOL.  XIV  ATT  1 
requirements  for  payload  bay  doors  open  operation 
2b.  Radio  wave  transponder  is  not  powered  during  launch  and  landing  and 

three  independent  inhibits  prevent  transmitter  operation  (see  attached) 


SAFETY  VERIFICATION  METHODS: 

1 .  System  test 

2.  Analysis  and  test 


STATUS : 


Open 

Approval 
Phase  1 
Phase  2 
Phase  3 


Payload  Organization 


STS 
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Inadvertant  turn-on  of  the  transmitter  is  precluded  by  three  independent 
inhibits.  There  are  three  relays  which  must  be  closed  before  the  transmitter 
operates  (see  figure  below).  The  relays  are  controlled  by  two  switches  of  the 
PSSP. 


Detailed  specification  of  the  transmitter  characteristics  is  ongoing. 
Currently  two  different  transmitter  designs  are  being  analyzed  as  indicated 
below. 


Design 

Operating 

Frequency 

Output  Power 
(Watt) 

Antenna  Gain 
(dB) 

Antenna  Size 
(cm) 

Baseline 

1-10  mHz  Swept 

1-10 

1-5 

50 

Alterna te 

5  GHz  Fixed 

1-10 

30 

15 
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PAYLOAD  HAZARD  REPORT 


NUMBER:  F5  ^ 


PAYLOAD: 

CIV 

PHASE : 

0/1 

SUBSYSTEM: 

Electrical 

Collision,  Fire, 
HAZARD  GROUP:  Corrosion 

DATE: 

Aug  87 

HAZARD  TITLE:  Battery  Explosion  or  Release  of  Corrosive  Materials  into 
Payload  Bay 


APPLICABLE  SAFETY  REQUIREMENTS:  NHB  1700.7,  Para.  209  HAZARD  CATEGORY: 

X  Catastrophic 
Critical 


DESCRIPTION  OF  HAZARD: 

Explosion/contamination/corrosion  of  payload  bay 


HAZARD  CAUSES: 

1.  Excessive  internal  pressure  due  to  high  temperatures 

2.  Under  or  reverse  voltage  in  cells 

3.  Overloads 

4.  Material  incompatibilities  or  failures 

5.  Environmental  loads 


HAZARD  CONTROLS: 

1.  Battery  pressure  re lief -thermal  analysis  will  be  performed  to  insure 
adequate  venting  under  all  thermal  conditions 

2.  Automatic  switch  off  in  undervoltage  condition 

3.  Internal  and  external  fusing 

4.  Selection  of  compatible  materials  or  adequate  protection 

5.  Compliance  with  structural  integrity  criteria 


SAFETY  VERIFICATION  METHODS: 
Analysis,  inspection,  test 


STATUS : 

Open 

Approval _ Payload  Organization _ STS 

Phase  1 _ 

Phase  2 _ 

Phase  3 
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